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SYHOPSIS 


A computational capability for the automated optimum 
elastic design of general framed structures is developed in 
the present work. The optimxmi design of cable stayed struc- 
tures, v/hich are becoming q[uite popular as long span bridge 
girders is studied using the developed programs. Ea,rlier studies 
in the analysis of this class of structures report that they 
exhibit geometrically nonlinear beha,viour under working load i 
due to the beam-column effect. The optimm design formulations 
are generalised to include this nonlinear behaviour^ I 

The objective function chosen is the total cost of materi 
used, which is to be minimised subject to constraints on stress^ 
and deformations in the structure. The bending elements of the: 
class of structures studied are assumed to be box sections of 
specified width and depth. The thicknesses of the wobs and 
flanges of these elements are the design variables. Constraints 
on the extreme fibre stress due to axial force and bending and i 
on the average shear stress in the webs are considered. These ; 
elements are supported at discrete points by high tensile I 

cables under an initial tension. The areas of the cables and 
the pretension forces are also considered as the design mriabls 
and constraints are imposed on the stresses in the cables due M 
initial tension and the stresses at working load. 
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The structural optimisation problem, y/hen cast as a non- 
linear programming problem, can be formulated in t?/o different 
ways: (1) Synthesis Approach and (2) Integra.ted Synthesis 
Approach. In the synthesis approach the geometric properties 
of the structure are considered as design variables, and cons- 
tieiints imposed on stresses and deformations. The structural 
analysis is performed explicitly to q.uantify the constraints 
by any of the well-known methods of structural analysis. The 
stiffness method of structural analysis is used in the present 
work. In the second approach, the eriuilibrium equations are : 
appended as equality constraints, considering the set of dis- ; 
placements as additional design variables. This eliminates the; 
solution of simultaneous equations overytime the objective 
function and constraints arc evaluated. Such a formulation is 
possible for discrete structures like trusses and frames using : 
the stiffness method of analysis and for continuous structures : 
analysed using the finite olemcnt displacement method. The syn- ; 
thesis approach leads to a mathematical programming problem witi 
nonlinear inequality constraints, and the integrated synthesis | 
approach to a problem with both equality and inequality cons- 
traints. The objective function in both the formulations turns ; 
out to be a linear function of the design variables. The non- | 
linear programming problem is solved using the Sequential Uncons 
trained Minimization Technique, using a quasi-Uov;ton algorithm 
for unconstrained minimization. 
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The nonlinear heam-col'umn effect in the frames is incorpc 
rated through the stiffness matrix of the elements* The secant 
stiffness matrix based on the conventional beam-column approach: 
is used in the present work. 

The computation of gradients required in the optimization 
procedure is generally performed by the finite difference tech- 
nique, This procedure, however, involves a considerable amount 
of unnecessary and repetitive computations. Efficient procedure 
are developed and adopted for gradient computation in the preset 
work. 

The two formulations mentioned above are applied to the 
minimum cost design of cable stayed structures in order to compa 
thoir efficiencies. Tbo synthesis approach has been found to be 
quite efficient and reliable in the present study. The intograt 
synthesis approach, on the other hand, suffers from the dis- 
advantage of proper scaling of the constraints to effect a smoot 
behaviour. Premature termination of the algoritlmi is experience 
for large problems and this can only be avoided by a proper 
choice of the scaling parameters. j 

The results of the example problems considered indicate ; 
that the minimum cost design exhibits the property of ’fully 1 
constrained design*. In other words, each geometrical design 
variable in the optimum design takes a value so that one of the | 
stress constraints of the corresponding member is active or tbe i 
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varialDle itself is at its lower bound, The optimum design, 
therefore, has at least as many bounded constmints as the 
number of geometric variables in the problem. The pretension 
forces in the cables alone are the free variables, and they 
govern the cost of the fully constrained design. Therefore, 
a third formulation is considered in the present work, in 
which the pretension forces in the cables are the only design 
variables, and at every design point, the fully constrained 
design is performed to find the cost. This formulation leads 
to an unconstrained minimization problem. The fully constrained 
design formulation is tho most efficient in the sense that it i 
is much faster than the synthesis ax)proach. ■ 

The effect, on the optimum design, of varia,tion of some 
of the input parameters, o.g. the ratio of cost of high tensile i 
cable to that of mild steel, minimum thickness requirements, 
depth of the girder and the allowable stresses in tho cable, I 
arc studied for both the linear elastic behaviour and the geo- 
metrically nonlinear behaviour and tho results discussed. | 

Some of the salient condusiono are: (l) The optimum I 
design of the structures considered is a fully constrained desigj 
in the sense defined earlier, and (2) for structures loaded j 
only at the nodes, inclusion of geometrical nonlinearity reduces; 
the minimum cost of the structure, while for structures loaded | 
along the members, an increase in the optimum cost is experiences 



The principal features of the present work are; 

(1) to compare the efficiencies of various optimal design 
fomiulations for ca,ble stayed structures, and (2) to emhed 
the geometrically nonlinear behaviour of framed structures 
in the optimum design framework and to study its effect on 
the optimum design of cable stayed structures vdiich have been 
reported to exhibit this behaviour. The work reveals that 
this behaviour docs not significantly influence the minimum 
cost design of this class of structures. 



CHAPTER ONE 


INTRODUCTION 


1.1 GENERAL REMARKS: 

Structural optimization may te defined as the ra,tioml 
evolution of a structural design which satisfies a set of 
functional requirements, and which, in terras of a predefined 
criterion, is the best of all such acceptable designs. Early 
attempts at the evolution of optimim designs led to the theorems 
of Maxvi?ell | lj and Michell ] 2j , and to considerable amount of 
analytical work in the recent years. Theorems providing the 
best choice of configuration and disposition of materials for 
structural elements and systems have been stated and proved. 

Ihe excellent surveys by Fasiutynski and Brandt j3| in 1963 
and by Sheu and Prager jAj in 1968 provide a complete historical 
development of the analytical methods of structural optimiza- 
tion. The developments in the area in the past seven years 
warrant another critical review. 

The advent of digital computers and their use in 
solving structural analysis and dosign problems led to the 
approach of Automated Optimum Structural Design, Developments 
in computer technology, providing the researchers with machines 
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of greater speed and storage capacity, loa-ve made this approach 
applicable to problems of varying degrees of complexity. Gon- 
seq_xiently, the methods of mathematical programming gained con- 
sidemble impetus and their use to solve problems of different 
degrees of complexity has, by nov;, become widespread. The 
latest innovation with on-line computing systems is the inter- 
active design in which the engineer is in constant communication 
with the machine during the design process. This development is 
yet in its infancy, and the current research in stimictural 
optimization is predominantly in the area of automated optimum 
design. 

A largo number of structural systems are idealised, 
analysed and designed as plane frames. Cable stayed structures 
form one such class of systems, which are quite popular among 
civil engineers. They aru mainly used in bridges, in the form 
of cable stayed girders, and in communication systems, as 
guyed towers. The main structural members are supported at 
intervals by high tensile wires under some initial tension. 

The cables being axially loaded members, the structural material 
is most efficiently used in these elements, and hence the 
material consumption of those structures is quite low as com- 
pared with structures of other types. 

The design of framed structures follows one of the 
following three design philosophies: a. Elastic Design, 
b. Plastic Design, and c. Limit State Design. Reinforced 
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concrete frames are designed by any one of the three philoso- 
phies, while steel structures are designed by either of the 
first two. Ihe study of optimum design of framed structures 
has received considerable attention of the researchers in the 
past, a,s is evident from the enormous amount of literature 
available today (5,6j. Curiously enough, attention was pri- 
marily directed towards the plastic design of frames by linear 
programming, and the elastic and plastic design of simple 
frames and their components by variational methods and conse- 

q_uent optimality criteria. Relatively loss attention has been 
paid to the elastic design of framed structures. I'he a.nalytical 
(variational) approaches are useful in arriving at qualitative 
conclusions and optimality criteria for different cases. They, 
however, hocomo extremely complicated and the solution seeking 
and interpretation of thu optimality criteria become difficult 
for any but the simplest of structures. The efforts have 
therefore been restricted to the design of simple beams, columns 
and portals with continuously or discretely varying sections 
subjected to simple loads. numerical methods, on the other 

hand, employ iterative or mathematical programming techniques, 
or a combination of both and they are used to design complex 
structures. Ibwcver, qualitative conclusions for design cannot 
be directly obtained by the numerical methods. They are popular 
because they are the only means of solving practieal problems of 
varied degrees of complexity. The analytical treatment of 
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simple structures, however, is necessary for the following 
reasons; It reveals faults in intuitive design criteria that 
have long heon tahen as self-evident (c.g. the simultaneous 
failure mode theory and the fully stressed design concept). Also, 
it provides a deeper insight into the optima, lity criteria, and 
era hies the designer to choose a near-optimum design as a 
starting point for the optimisation using thu numerical tcch- 
niduos . 

The plastic design problem, under certain simplifying 
assumptions, can bo ca:St as a linear programming problem, which 
can bo solved using standardised procedures. The major draw- 
back of the general plastic design formulation is tha.t the 
serviceability constraints under working load arc not considered. 
The inclusion of these constraints would make the formulation a 
nonlinear programming problem. The elastic design invariably 
loads to a nonlinear programming problem which has no specific 
standardised solution procedure. Most of the effort in the 
elastic design has been directed towards developing efficient 
procedures to solve the resulting nonlinear programming problem. 
The complexities of nonlinear programming problem are probably 
the deterrants which divert the rosearchors away from elastic 
design, although it is one of the most challenging fields of 
optiiml structural design. 
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The present work addresses itself to the automated 
optimum elastic design of cable stayed bridge girders. Three 
formulations leading to nonlinear programming problems have 
been studied and their relative efficiency discussed. The 
effects of inclusion of the geometrically nonlinear behaviour 
of the structure, and of the variation of some of the predeter- 
mined design parameters, on the optimum design, are studied. 

1.2 THE IRESENl TORK: 

The computational capability for the automated optimum 
design of planar framed steel structures under linear elastic 
behaviour is developed in the present work. The objective 
function chosen is the total cost of materials used, which is 
to bo minimised subject to constraints on stresses and deforma- 
tions in the structure. The design variables chosen arc the 
cross sectional properties of the members. 

The structural optimization problem, when cast as a 
nonlinear programming problem, can be formulated in two different 
ways: one in which the geometric properties of the structure are 
considered as design variables, and constraints imposed on 
stresses and do fo mat ions. The analysis is porfonnod explicitly 
by any of the well-known methods of structural analysis. The 
stiffness method of analysis is adopted for this formulation 
in the present work. This formulation, henceforth, is termed 
the ’Synthesis Approach', the details of v/hich are presented in 
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Section 2.2. Alternatively, the equilihrinm equations are appends 
as equality constraints, considering the set of displacements as 
additional design variables. This eliminates the solution of 
simultaneous equations everytime the objective function and 
constrt\ints are evaluated. Such a formulation is possible for 
discrete structures like trusses and frames using the stiffness 
method of analysis, and for continuous structures idealised and 
analysed using the finite element displacement method. This 
formulation, henceforth teraiod the ’Integrated Synthesis 
Approach', is described in Section 2.3. The synthesis approach 
leads to a mathema.tical programming problem with nonlinear 
inequality constraints and the integrated synthesis approach 
to a problem with both nonlinear equality and inequality cons- 
traints. The objective function in both the formulations turns 
out to be a linear function of the design variables. 

The following are the general features of the computer 
programs developed for the two formulations mentioned above. 

They are independent of the structural configuration, in the 
sense that all the relevant data required for the optimization 
is automatically developed from the geometrical and other data 
of the structure, provided as input. The objective function, 
constraints and the gra,dieuts with respect to the design 
variables are evaluated implicitly at any design point. They 
also admit the possibility of different members being made of 
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different materials with different unit costs, and also of 
different shapes. Given the type of member, the program 
automatically chooses the design variables and the stress 
constraints applicable to that type of member. Such a capability 
has also been incoirporated earlier in the works of Kavlie and 
Moe j7| for the optimum design of the component fmmes of 
ship structures. 

I'hc present work aims to study the optimum design of 
planar framed steel structures supported at discrete points 
by high tensile cables under initial tension, vrith particular 
rofcrenco to cable stayed bridge girdors. The bending elements 
in these structures aro subjected to a considerable amount of 
axial force and earlier works in the analysis of these struc- 
tures report that they exhibit geometrically nonlinear behaviour 
due to the beam-colTjmn effect oven under working loads. Ihe 
above formulations are therefore generalised to incorporate 
this nonlinear effect with little modification to the computer 
programs developed. The relative merits of the two formula- 
tions arc studied through the example problems. 

fhe example problems solved during the course of this 
work show that the optimm design for this class of structures 
is a 'Fully Constrained Design'. In other words, each geometri- 
cal design variable in the optimum design takes a value so that 
one of the stress constraints of the corresponding member of 
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the structure is active or the variahle itself is at its 
lower bound. Ihe optimum design, therefore, has at least as 
many bounded constraints as the number of geometrical design 
variables in the problem. The pretension forces in the cables 
are the free variables, which control the behaviour of the 
structure. Therefore, a third formulation is considered in 
the present work, in which the pretension forces in the cables 
are the only design variables, and at every design point the 
'fully constrained design* is performed. This formulation 
leads to an unconstrained minimization -problem presented in 
Section 2.4. The computer program developed for this formula- 
tion also preserves the generality and the versatility of the 
earlier two. 

General approaches to the analysis of geometrically 
nonlinear plane frames and the modifications required in the 
computer programs to incorporate the nonlinear effects are 
discussed in Chapter 3. The nonlinearity in the present -work 
is included in the stiffness matrix based on the conventional 
beam - column theory. 

The comiputation of gradients of the objective function 
and constraints for such general formulations is usually per- 
formed by the finite difference scheme. However, the finite 
difference procedure Involves a considerable amount of unneces- 
sary and repetitive computations. Alternative approaches to 
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the problem of rearalysis of modified structures, and the 
efficient procedures developed and adopted for gradient 
computation in the present vrork are detailed in Chapter 4. 

The formulations are studied using the examples of a 
tied cantilever and two cable stayed girders subjected to 
single load condition. Results are obtained for both the 
linear elastic behaviour and for the geometrically nonlinear 
behaviour of the structure. Also, some of the input design 
parameters are varied and the effect of the variation of 
these parameters on the optimum design are studied. The results 
of these studies are presented in Chapter 5 and the discussion 
of the results, the conclusions drawn and the suggestions for 
extension of this work are presented in Chapter 6. 

1.5 BEYIEW OR CURRENT IITEMTURE II T,HS RELATED TOPICS: 

'One of the first attempts at the optimum clastic design 
of framed structures was made by Brown and Ang ] 8] , who applied 
the gradient projection method to the design of frames. A 
continuous relation between the area and moment of inertia of 
wide flange beam members was determined, and each member was 
defined by a single design variable, its moment of inertia. 

They have shown that the minimum weight design under multiple 
loading conditions, is not a fully stressed dosign, Reinschmidt 
ct al joi have shown that the general iterative methods need 
not necessarily converge and even if they do, it need not 
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nccGssarily bo to tho minimtiin weight design. The successive 
linearisation approach was proposed by them for the design of 
frames. Rorastad and Yfeng | loj solved tho nonlinear programming 
problem by successive linearisation and application of linear 
programming. 

Kavlie and Moo jllj and Moses and Onoda ll2j hawc consi- 
dered tho problem of minimum woight design of elastic grillages. 
The problem is a nonconvex programming problem, which is solved 
by Ka,vlie and Moe by the penalty function approach. 

Majid and Elliott j 13] considered the optimum design of 
a gable frame subjected to deflection constraints, and solved the 
problem by a 'dynamic search method’. They have presented 
design charts for the optimum design of gable frames of arbi- 
trary span and height subjected to the deflection specifications 
of BS 449. The formulation does not consider stress constraints. 
Majid and Anderson jl4j formulated the x->noblem of optimum design 
of general indeterminate structures using the force method of 
analysis. The constraints derived had separable variables, and 
they vi-ere linearised piecewise, and solved using the simplex 
method, 

ICavlic and Moe j7| have employed the seq.ucntial uncons- 
tmined minimisation technique for the optimum design of framed 
ship structures subjected to different loading conditions, - ■' 
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Tbeir formula -tion is very general and ha.s been quite successful!; 
applied to the design of different structures. One formulation ; 
in the present work, namely the synthesis approach, has a cer- 
tain amount of similarity with this v/orlc in its generality and 
the method of solution. 

A different approach has also been tried by some authors. 
The conventional mathematical programming approaches to the 
optimal structural design suffer from the disadvantage that a 
number of analysis - design cycles is necessary to arrive at 
the optimum design. A . possible way of reducing the computationa: 
effort is to combine the mathematical programming methods with ^ 
some iterative techniques. Wright and Peng jl5l considered 
such an approach and termed it the 'multimode approach’. In 
this scheme, the stress ratio iteration method ll6j, the ba.sic 
iteiation method j9l and the gradient projection method are ' 

used sequentially. The use of three different techniques sequen-: 
tially results in the exploitation of the relative merits of ' 
each of these techniques to form an efficient composite method. 
Also, a detailed derivation of the behaviour gradient matrices 
is presented by them; these are used to compute the rates of 
change of the constraints with respect to the design variables, 
lavies and Wang j 17l used a combination of iterative fully 
stressed design and linear approximation tochniquo to obtain 
the optimum design of railv/ay underframos, idealised as space 


frames, 



12 


The very challoiiging problom of optimization using 
discrete available sections has been considered by Celia and 
Logcher | 18j . An overview of the algorithm for discrete 
optimization and its application to some practical design 
examples are presented by Celia and Soosaar !l9j. 

Gable stayed structures, especially for medium span 
bridges, have been quite popular in Europe for quite some time, 
and recently Leonhardt and Zellner j 20| iB ve shov/n tha.t cable- 
stayed girder bridges are superior, in terms of economy and : 

structural efficiency, to other types of bridges for all spans 
over 200 metres. The analysis of these bridges was considered 
by Stafford Staith j2lj who used a mixed force-displacement 
approach, which he extended later {22{ for the analysis of | 

double plane cable stayed structures. Troitsky and lazar j23j 
used a flexibility approach while Lazar [ 24j used a stiffness ; 
approach. It was observed that the designs exhibited geometri- 
cally nonlinear behaviour, which was considered by Lazar {24| 
who used a Nowton-Raphson type solution procedure, Tang j 25l 
used a transfer ma.trix technique for the analysis, including 
the nonlinear terms in the form of imaginary external loads. 
Kajita and Cheung | 26| performed the linear static and dynamic 
analysis of the structure by the finite element method, consider- 
ing the deck as a plate and the tower as a beam. Morris j 2?! 
performed the dynamic analysis of the planar structure both under 
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linear and nonlinear letoavionr, and found that the effects of 
the nonlinearities 'were marginal. 

An approach to the efficient design of these struc- 
tures is presented by Lazar et al ] 2S| v/ho present a kind of 
load balancing approach in which the cable pretensions (or 
post- tensions as they call it) are chosen in such a ■v'/a,y that 
the maximum bending moment along the girder is minimised by 
an iterative technique using the flexibility method. It may 
be mentioned that the choice of minimum design bending moment 
along the girder as the optimality or efficiency criterion is 
rather dubious. Since the structure is statically indeter- ; 

minatc, the properties of the sections would have to be deter- 
mined before the analysis, Ihe formulations considored in the 
present work include the tensioning force as an additional 
design variable, thereby controlling the behaviour of the 
structuro. Tang |29l highlights some of the prohloms encountered 
in the practical design of these structu.res, Ibwovor, no effort 
seems to have gone to seek the optimal design of this class of 


structures. 



CHAPTER TW) 


THE PROHLBM POEHUIiATIOH 


2.1 IHTRDDUGTIOH: 

The general problem considered in the present work is 
the optim-uir design of planar steel frames with linear and 
nonlinear elastic behaviour under single loading condition, 
with special reference to cable stayed planar frames. Both 
joint and member loads are considered. Provision is made 
for different types of member cross section. 

The different types of members considered in the 
programs are: a. Box beam, b. I-Section, c. latticed section of 
triangnlar or square shape, and d. High tensile cable. Any 
other structural section can be included with little effort. 

For each type of member there are some preassigned parameters 
which remain constant during the process of optimization and 
some variable parameters which form part of the design vector. 
The members are assumed to be prismatic. Stress constraints are 
imposed on each member and lower bounds on each variable. A list 
of preassigned parameters, design variables, expressions for the 
cost of the member, and the stress constraints considered for 
each member type included in the computer program is given 
in Tabic 2.1. 
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The objective function chosen is the cost of materials. 
The total cost is computed as the sum of the costs of consti- 
tuent members. The cost of a member is simply the volume 
times the cost per unit volume of material for that member. 

The inclusion of fabrication and erection costs of each mem- 
ber, and the costs of additional structural components like 
the web stiffeners of beams and bracings of latticed members 
can be admitted in an indirect manner. This is sijnply achieved 
by increasing the unit cost of the member by a certain per- 
centa,ge . 

from Table 2.1 it may be observed that the objective 
function is a linear function of the design variables for the 
different cross sections considered. The depth of the box 
beam or I-soction is considered as a proassigned parameter 
because it is quite often determined from nonstructural consi- 
derations. Purthermore, for large span bridge sections, 
shapes are determined from aerodynamic stability considera,tions , 
and this results in fixed values of the width of flange and 
depth of the idealised cross section. 

2.2 TIffi SYNTIIESIS APPROACH 

2.2.1 Amixaia: 

This formulation uses the stiffness method of analysis. 
The stiffness matrix of the structure is generated at every 



18 


desigri point during the search procedure, and the equilihrium 
equations are solved to obtain the joint displacements, Ihe 
nonlinear effects are embedded in the stiffness matx'ix formu- 
lation as described in Chapter 3. Once the joint displace- 
ments are known the maximum forces in the members are evaluated 
and then the maximum stresses are computed. The relevant 
equations are given below: 

'ihe equilibrium equa.tions: 

K u = J (2.1) 

in v/hich, 

K is the stiffness matrix of the structure, 
u is the nodal displacement vector, and 
p is the nodal force vector. 

Member end forces for member i: 

f. = S. T. 6. + fi. (2.2) 

in which, 

f is the vector of member end forces in tho olemicnt 
coordinate system, 

5 is the stiffness matrix of the member in the 
element coordinate system, 

T is the coordinate transformation matrix, 

6 is the vector of joint displacements of the two 
ends of the member in the global coordinato system, 
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ff is the vector of fixed end forces in the member 
in the 'element coordinate system, 
and the subscript i denotes that these quantities are related to 
member i. 

The maximum bending moment, if any, along the length 
of the member is computed from statics. The design bending 
moment, shear force and axial force for the member are then 
chosen and the constraints evaluated. 


2.2.2 C Qnstra.ints : 


The maximum fibre stress in combined bending and axial 
force is given by, 


(T. 

X 


Q- 

A. 

1 


+ 


M. 

1 

2i 


(2.3) 


in which, 

ar is the maximum fibre stress, 

Q is the axial force, 

A is the area of cross section, 

M is the design bending moment, 

Z is the section modulus of the member, and 
the subscript i denotes that these quantities are related 
to member i. 


The average shear stress in the web is given by, 

Y. 

q. = 

1 Aw^ 


(2.4) 



120 


in which, 

q_ is the average shear stress, 

V is the design shear force, and 
Aw is the area of the wet. 

The stress constraints are written in the form 

1 20- ^2-5) 

and 1 - (ili/qaii,i) >. 0 (2.6) 

O'-,, . and Cl T T . are the allovra.ble maximum fibre 

d JLX ,1 3.11 5 1 

stress and shear stress respectively for member i. Normaalization 
of the constraints, as indicated in Eq_uations 2,5 and 2.6, faci- 
litates a smooth behaviour of the optimization problem. The 
constraints would always lie between 0 and 1 , and the value of 
nearly 1 for a constraint shows that the corresponding stress 
is nearly zero in the member. 

For high tensile cables, different sets of constraints 
are used. The limitation on the pretension stress is given by 

P. 

1 - ( 3^ ) / (0.6'PP >. 0 (2.7) 

in v-^hich, 

P is the pretension force in the cable, and 

rr is the yield stress or the proof stress of the 

y c 

material of the cable. 
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The limitation on the- fina.l stress in the oahlo is given 




1 - (^) / (0.8 o: ) > 0 


^A. 

1 


yc^ 


( 2 . 8 ) 


Displacement constraints, if any, are appended to the 
above constraint set, 

2.2.3 T he Opt i mization Proble ra; 

The design vector in this formulation consists of the 
member dimensions and pretension forces. The formulation 
results in the follovang mathematical programming problem: 


Dind X a. d to minimise f(d) 
Subject to 


mi 
£ C. 
D=1 


D 


g^('i) > 0, j = 1,. , DC 


^i 

u?- 


> q, 


1 ,..., im 


(2.9) 


< Ui. < u^!: for k 6 M 


"k -• “"k - nc 

in v/hich, ' , 

X is the : .general design variable vector of dimension D, 
d is the subset of x containing elements corresponding 
to member dimensions and pretension forces of 
dimension NDV. In this formulation x =: d, and 
M = NDV. 

C. is the cost of the j-th member, 

J 

It.i is the number of members, 
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g. is the j-th iDehaviour coiistxaint , 

J 

NC is tho nmnbGr of behaviour constraints, 

is the lower bound on the i-th design variable , 

1 H 

and u^ arc respectively the upper and lower 
bounds on the k-th displa coiacnt , and 

M is the number of degrees of freedom of the structure. 

A flow diagram of the procedure for evalmtion of the 
objective function and constraints in this formulation is given 
in fig. 2.2. 

2.3 im lUTEGEAIED SYNltlESlS ilPPEDAGh: 

Ihe synthesis approach recLuires that the structure be 
analysed by solving the eq.uilibri'um equations at every design 
point during the search procedure. This consumes a considerable 
portion of the computational effort. The 01113/- waj to avoid the 
solution of simultaneous equations so many times is to consider 
the displacements of joints as additional design variables and 
to treat the equilibrium equations as equality constraints. 

Such a formulation, of course, increases the mrnibor of variables 
and constraints and introduces the troublesome nonlinear equality 
constraints to the optimization problem. The purpose behind 
the study of this formulation is to compare the efficiencies 
of the two formulations. It may not be out of place to mention 
here that a similar formulation for the optimum, design of truss 
structures was made by Schmit and Fox l30,3ll. 
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The design variable vector, thorofore, consists of two 
parts - the joint displacements, u, and member dimensions, d. 


. * . X 

2.3.1 Constraints : 



( 2 . 10 ) 


The inequality constraints on stresses arc computed from 
the displacements in the same manner as detailod in Section 2.2.2. 
The equality constraints on equilibrium arc written as, 

K u - p = 0 (2.11) 

The elements of the stiffness matrix, K arc linear 
functions of the dimension variables d. The constmints given 
by Eq. 2.11 arc, therefore, nonlinear equality constraints as 
products of dimension variables and displacement variables 
appca3: in those equations. To effect sraooth behaviour of the 
optimization procedure, each of the equations is noiraalised 
with respect to a prcsuloctcd force value. Tot this value, 
bo represented by EHA.X^ for the k-th oq.iiilibriurn equation. The 
k-th normalised equality constraint therefore becomes, 

Ik "^r " ^k 1 ^ (2.12) 


The preselected normalising force and moment for the 
structure arc fed in the computer program and the value 
is picked from those two depending upon whether the k-th equation 
expresses the force equilibrium or the moment equilibriun. 
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2.3.2 The Optimization Probl om: 


This prohlem thus results in the following nonlinoar 


programming prohloin: 


Pind J 


4 


^ IM 
to minimisG f(d) = 2 

D=1 


Subgoct to gj(x) 2 0? j = 1,..., KG 


l-[^(x) = 0, k = 1 , M 



2 

d. - 

1 


1 ^ 

1 

^k 

i. 

^k 


.u 


(2.13) 


In this formulation the total number of design variablos 
N = mi + M. 


A flow diagram of the procedure for evaluation of the 
objoctivo function and constraints in this fomiulation is 
given in Pig. 2.3. 

It may be mentioned hero that the intormodiato designs 
obtained during the optimization procedure in this formulation 
will not be physically meaningful, since they v;ould not satisfy 
the equilibrium equations. 


2.4 Tl-iS PUIIY COISTRAINEI nesIGU PORMTJMTION: 


The optimal design results of the cable stayed planar 
frames studied by the above two formulations indicate that the 
optimum design of this class of structures is a 'fully constrained 
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design’, that is, every design variable defining the cross 
sectio’nal property of the elements takes a value at the opti- 
mum design such that it either satisfies a stress constraint 
or is hounded hy the side constraint. In a box beam or I-Scction, 
the thickness of the web is governed either by the minimum 
thickness constraint or by the shear stress constraint. The 
thickness of the flange, likewise, is governed either by the 
extreme fibre stress constraint or by the minimum thickness 
constraint. In the high tensile cable, the area is governed 
cither by the ollov/ablc pretension stress or by the allowable 
final (or v/orking) stress. Thus the total number of active 
constraints at the optimum is at least equal to the nuiabor of 
dimension variables defining the cross sections of the members. 

The pretension forces in the cables, hov/ovor, arc the only free 
parameters which influence the cost once the constraints arc 
satisfied. 

The iterative procedure of a 'fully constrained design' 
in which members are resized by the satisfaction of the cons- 
traints suggests itself as an efficient design procedure. The 
pretension forces can be taken as design variables in an 
optimization formulation to minimise the cost of materials. The 
design would be performed as a ’fully constrained design' at 
every design point. The behaviour constraints and the side 
constraints are implicitly satisfied in this approach and 
hence the resulting optimization problem is an unconstrained 
minimization problem. 
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Tbe relevant eq_'aations used for "the design of components 
by the fully constrained design procedure are given below: 

a. Box beam and I-Section: 

file area of the web is given by A, = 2 Dt, (2.14) 

V/ w 

I'o satisfy the shear stress constraint, the 
required v/eb thicloness is , 


t 


Y 


(2.15) 


W1 

Iherefore, the thickness of the web is chosen as 

t = MAX (t ,, t ( 2 . 1 6 ) 

where t is the mininiurfl allowable thickness of tho web. 


wmin 

The area of cross section is given by, 


A = 2(Bt„ + Dt ) 

^ f vr 


.(2.17) 


and the section modulus is, 

Z = D (Bt^ + Dty3) (2.18) 

'I'ho combined bending and axial stress is given by Equation 2.3. 
Since t^^ is known, the required thickness of tho flange t^^ 
is obtained by satisfying this stress constraint, viz. , 



Sl 


2(Bt„. + Dt ) 
^ f1 w^ 


4 * 


M 


D(Bt^^ + Dty3) 


(2.19) 
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'I'he thickness of the flange is then choson as 

= MAI (tj^, .(2.20) 

where is the minimum allowahle thickness of the flange, 

h . H igh Sonsile Ca.blG : 

Iho area of the cable is obtained from, 


= MAX ( 


.A 


0.6 <r ’ O.S(r ’ '^cmiii 
ye ye 


,, ) ( 2 . 21 ) 


whore minimum allowable area of the cable, 


c. .Iria, iig;ula..XL.. La tt i co ; 

I'hc required area of each log of the lattice, A^.|j is 
found from the satisfaction of the corebined axial and bending 


stress, 


’^11 ^ ^ 1 IL 1 (2.22) 

in v/hich LI is the lateral distance between the legs of the la 
lattice. Ihc area of the leg is chosen as, ’ 

-1 = (Aj^, ) (2.23) 

where is the minimum allowable area of the log. 

d. battiooi : 

I'ho required area of the log in this ease is found from 


(T 


all 


4A. 


'L1 


•f 


M 


2Aj^ ^ . IL 


(2.24) 
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and the area of "tho log is choson using Equation 2.23. 

2 . 4- . 2 The Ontimiaa tion Problem ; 

The fully constrained design leading to the minimum 
cost design formulation, therefore, results in an unconstrained 
minimization prohleni, with the iterative fully constrained 
design performed at every stage. Thus this fomulation can 
he considered a combination of the mathematical programming and 
optimality criterion techniques. A detailed flow chart of the 
optimization process is given in Fig. 2.4- and a flow chart for 
the design is given in Fig. 2.5. 

This formulation cannot directly handle displacement 
constraints. This is because the optimum design of a structure 
with active displacement constraints is, in general, not fully 
constrained in the sense defined earlier. 
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Generate element stiffness matrix and transfounation 
matrix for the memher. 

Find member forces in element coordinates 
Find maxim-uon B.M. , shear and axial force in 
the member 



Fig. 2.2: Flow Chart for svaluation of Objective Function 
and Gonstraints-Synthesis Formulation. 
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Pig. 2.3: Plow Chart for Evaluation of Objective Punction 
and Constraints - Integrated Synthesis Approach 
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Blow Chart for Design and Evaluation of 
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CHAPTER THREE 


NOKIilNEAR ANALYSIS 


3.1 I'NTRODUGTION: 

The norilinear beha viour of structural systems is broadly 
classified into three categories: (i) Geometric Nonlinearity, 
which arises from the nonlinear terms in the kinematic equations, 
(ii) Material Nonlinearity, which arises from the nonlinearity 
in the constitutive equations, and (iii) Combined Material and 
Geometric Nonlinearity, 

The structures considered in the present study exhibit 
only geometric nonlinearity, since at the working load the 
raa,terial of the structure is assumed to be linearly elastic. 

These structures are supported by high tensile cables with a 
certain amount of initial tension and the flexural eJ.ements 
are subjected to a significant amount of axial force. The 
presence of axial force reduces the flexural stiffness of these 
elements and increases the fixed end moments. This interaction 
between the axial load and bending is commonly knov/n as the 
beam-column effect and is the most important reason for the 
nonlinear behaviour of these structures. A certain contribution 
to the nonlinearity is also made by the large deformation of 
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structures, in which case the equilibrium equations should be 
written in the deformed configuration of the structure. This 
effect, however, is very small and the deformations of the 
structures considered in the present work do not wa.rrant the 
inclusion of this effect. The high tensile cables sag due to 
self weight and their behaviour is also nonlinearly elastic l20|. 
This nonlinearity is significant only for cables of very large 
spans and when the stress level in the cables is low (Figure 6 
of Ref. 20). The optimum design in the problems considered in 
the present v/ork results in the pretension and final stress 
levels to be relatively high and this justifies that the nonli- 
nearity due to cable sag be neglected. 

A brief review of the techniques of nonlinear analysis 
of structures is presented in this chapter. The choice of the 
'Secant Stiffness I/Iatrix' for the beam elements used in the 
present work is discussed, and the modifications to the computer 
programs necessary to incorporate the nonlinear analysis are 
also outlined. 

3.2 FOKEIilGAL TECHNIQUES FOR NOiOilHEAR STRUCTURAL ANALYSIS:. 

The techniques of nonlinear analysis followed the methods 
of linear analysis and naturally nonlinearities wore attempted as 
corrections to the linear solutions. The methods developed from 
these were generally iterative, incremental or a combination of 
both. Fig. 3.1 shows the Newton-Raphson or the 'Ta-ngcnt Stiffness' 



FIG 3'1 TANGENT STIFFNESS 
APF’l-iOACH 


FIG.32 SECANT STIFFNE 
APF^ROACH 


Iteration 


FIG. 3-3 PURELY INCREMENTAL 
APPROACH' 


FlG.3-i INGREMENTAL-CUM- 
ITERATIVE AF'PROACH 


METH0D50F NONLINE AR STRUCTlIKAl ANAIY‘,IS 
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approach to a typical single degroo of freedom problem in struc- 
tura.l analysis. Fig. 3.2 shows the direct or 'Secant Stiffness' 
approach to the same problem. Both are iterative methods of 
analysis. The tangent stiffness approach is more efficient, 
in terms of rate of convergence, than the secant stiffness 
approach. Pig. 3.3 shows the purely incremental technique of 
solving the same problem. In this approach the error in the 
solution increases with the increase in the load. This is 
corrected by incorporating the tangent stiffness iterations 
after every few increments. Such an incremental-cum-iterative 
approach, shown in Pig. 3 . 4 ? is the most efficient and accurate 
technique for studying the nonlinear load-deformation history 
of a structure. 

Another t 3 rpe of corrective approach to tho linearised 
equations is the 'initial stress' or 'geometric stiffness' 
matrix. This arises due to the effects of the initial stress 
during an increment performing work on the incremental changes 
in the nodal displacements of the element. This approach is 
entirely satisfactory for the classical stability analysis of 
structures with negligible prebuckling deformations. However, 
it is insufficient for the cadculation of general geometrically 
nonlinear behaviour i32l. 

Other efficient methods for tracing the complete load- 
deflection history of highly nonlinear structures are based 
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on the parametric differentiation method. Self correcting | 531 
and numerical integration methods to solve the resulting initial 
value problem have been used by various investigators. 

Analysis of nonlinear structures can also be performed 
by the direct minimization of the potential energy functional 
134,35]. This approach, however-, suffers from the disadvantage 
that proper scaling of variables in the* analysis problem is 
req.uircd to ensure fast and accurate convergence. Computationally, 
this method is not competitive with the methods detailed 
above. 

All excellent comprehensive assessment of the solution 
techniguos applicable to static nonlinear structural anjilysis 
is presented by Tillerson, Stricklin and Haisler j 36j . 

3.3 NOlffillCEAR ANALYSIS APPROACH IN IHE PRESENT \OPJi; 

In the search procedures for optimization, a soguence 
of designs is generated such that the obyectivo function 
decreasos monotonically. Lot d^ and be two consecutive 

design points. At point i the cq_uilibrium equations are 
given by, 

[K(di, u^)] u. = 5^ (3.1) 

?/hich, lot us say, have boon solved to find u^. At point 
(i+1) the equilibrium equations are given by 

^i+1 = Pi+1 




(3.2) 
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The load vector also charges from point to point because the 
fixed end moments also depend on the axial load in the non- 
linear analysis. 

Let u^_^^ = + Ziu (3.3) 

lor the tangent stiffness approach we require the approximate 
estimate of the displacements and the unbalanced forces under 
these estimated displacements, let the approximate estimate 
of the displacements be u^- . The unbalanced forces at point 
(i + 1) under these displacements arc given by, 

= Pq+i - ^i <^5.4) 

The fengeiit stiffness equations are 

[TK(d^_i_^, u^)] ,Au = AP (3.5) 

in which TK is the tangent stiffness matrix at point i+1 evaluated 
using u^ as the estimated displacements, au is found by solv- 
ing the Equations 3.5, and thus is obtained. Thus in this 

approach it is necessary to evaluate two stiffness matrices, 

K to find the unbalanced forces and TK to find the corrective 
displacements. 


In the direct approach we assume that the design points 
i and i+1 are sufficiently close, so that the displacements 
u. and u. would also not be very different. In such a case 

1 1t1 

we may permit the approximation that 






(3.6) 
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Equation 3 . 2 may "be then approximated to 

Sp] = 5.^^ (3.7) 

The stiffness matrix is computed and Equation 3,7 solved to 
obtain directljr. This approach requires the assembly and 

inversion of one matrix, and is sufficiently accurate if the 
nonlinearities are small and the points are quite close. It 
is known that the structures considered in the present work 
exhibit only a small amount of nonlinearity. It is, therefore, 
decided that the direct formulation, as given by Equation 3.7 
would bo sufficiently accurate for embedment in the optimiza- 
tion procedure, furthermore, if the analysis is performed 
iteratively at the beginning of every linear search in the 
optimization procedure, any error in the displacements accumu- 
lated over the previous linear search would be Gllniin£i ted. 

The iterative equations are given by, 

[K(d, f = 5(7+1) 

in which the superscript in parentheses denotes the iteration 
count . 

The beam clement may be modollod either by the finite 
element approach | 37 - 45 j or by the stiffness approach using 
the conventional beam-column theory j46-54j. The finite 
element approach using the generalised polynomial requires 
nimierical integration to be performed to obtain the nonlinear 
stiffness matrix. The nonlinear stiffness matrix based on the 
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Toeam-coluinn theory is the easiest to program, and requires 
very little extra, computational effort over the linear stiffness 
matrix. In fact the linear stiffness rntrix is a particular 
case of the nonlinear stiffness matrix. Ihe best choice of the 
stiffness matrix is, therefore, the secant stiffness matrix 
based on the conventional beam-column theory. The stiffness 
matrix is given in Appendix A. 

It may be mentioned here that Equation 3.2 is directly 
used in the constraint set of the integra,ted synthesis approach, 
because in this formulation both d and u are the design vari- 
ables. Thus the nonlinear eqmlity constraints do not involve 
any approximation wha.tsoever. 



CHA.PTER POUR 


THE I'fflTHOD OP OPTIMIZilTION 

4.1 GI3DICE OP THE METHDU: 

Phe formulations of the optimal design prohlem considered 
in the present work lead to nonlinear programming problems with 
linear objective function and linear and nonlinear constraints. 

One of the objectives of the present work is to compare the 
efficiencies of the first two formulations described in ^ 

Chapter 2. This comparison should be done with respect to one 
optimization program. The three general classes of widely used 
nonlinear programming methods are: a) The Gradient Projection 
Method of Rosen |55| modified subsequently by Goldfarb (56l. 

This works well with linear constraints. However, the efficiency | 
is considerably reduced in the case of nonlinear constraints ; 

and nonlinear eq.nality constraints would be extremol 3 ^ difficult 
to handle, h) The Feasible Direction Lie thod of Zoutendijk \5l\ \ 

would result in a similar difficulty. The usual problem of j 

zigzagging would be more acute in the presence of nonlinear j 

equality constraints, c) The Barrier or Penalty function Methods | 

I 

arc more reliable, albeit slower. The main a,dvantago of using 
the penalty function methods is that the constrained minimize- | 
tion problem is transformed to a sequence of unconstrained 
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minimization problems which can he solved with relative ease. 
These methods are ofaite reliable and the constraint repulsion 
character of the methods makes the optimization process smooth. 
The main drawback of the method is that the unconstrained 
minimization problem has to be solved a nimiher of times and 
therefore the computational effort is more. Additional compu- 
ter time, however, is well worth the reliability the method 
generally offers. Another shortcoming of these methods is 
the severe ill-conditioning of the unconstrained minimization 
problem as the value of the penalty parameter tends to its 
limiting value (infinity or zero). Powell | 58j ha,s suggested 
a penalty function that avoids the ill-conditioning in eq.uality 
constrained minimization problems, and this has been modified 
and generalised by Osborne and Ryan |59| for both equality and 
inequality constraints. However, these proposals require 
extensive numerical experimentation and modifica.tion boforo 
they can be used for complex practica.1 problems. The penalty 
function suggested by Fiacco and McGomick j 6oj is chosen for 
the present work since it has been studied quite extensively 
and found to be reliable. The penalty function is of the form 

r) = f(x) + r 2 — + -^ 2 1^ (x) (4.1) 

D=1 g^Cx) 

in which, 

^ is the penalty function, 
r is the penalty parameter, and 

HOC is the total number of inequality constre.ints in the 
problem. 
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This is a.n interior penalty function approach and therefore 
the value of r reduces for successive unconstrained minimiza- 
tions. 

The subroutine VA10A developed by flotcher [ 61 [ is used 
for the unconstrained minimization. This program makes use 
of a q,uasi-Newton algorithm and follows the line of the 
Davidon-Fletcher-Powell algorithm. The Hessian matrix is 
approximated and updated during successive linear searches.. 

(In the DIP algorithm the inverse of the Hessian is a,pproxima,tcd 
and updated). The program uses a single quadratic fit for 
linear search. Some studios by Fletcher \& 2 \ indicate that 
quasi-Hewton methods are most efficient if the linear search 
is carried out v/ith minimum number of function evaluations. It 
is likely that such an approach may not be as robust on ill- 
conditioned problems such as those that arise when using penalty 
functions. This aspect is discussed in Chapter 6. 

4.2 CRADIENT COMPUTATION: 

The unconstrained minimization procedure used in the 
present work is a gradient method and requires the gradient 
of the penalty function to be evaluated at the start of every 
linear search. The gradient of the penalty function is obtained 
by differentiating Equation 4.1 with respect to the design 
variables , 



47 




+ 



M 

2 

T £=1 




i 


( 4 . 2 ) 


Thus the gradient computation requires the evaluation of the 
rates of change of the objective function and the constix'vints 
with respect to the design variables. These are usua,lly 
computed for such general optimization formn.lations by the 
finite difference scheme. The subroutine VA10A requires only 
about 3 to 4 penalty function evaluations per linear search, 
and the total number of linear searches is increased because 
the linear search is not very accura.to. Thus the use of this 
program for the optimization problem with 14 design variables, 
combined with a foiward or backward difference schcrac vrould 
require about N+4 penalty function evaluations per linear search, 
out of which N+1 would be required for the gradient computation 
alone. Thus any alternative procedure which leads to a faster 
coraputation of the gradient would improve the officioncy of 
the optimization program substantially. 


A major portion of the computational effort in the 
evaluation of the objective function and the constraints at a 
design point is spent in computing the stiffness matrices of 
the elements, and in the solution of the equilibrium equations. 
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The use of the finite difference scheme results in a major 
portion of this effort heing repoatod several times; this repe- 
tition, if avoided, can result in a considcrahle saving in the 
computer time.. Alternative procedures for gradient evaluation 
are developed and implemented for the synthesis and the inte- 
grated synthesis approaches. The procedures are hased essen- 
tially on the technirj^ues of roanalysis of modified structures, 
which arc discussed in the next section. The procedure for 
gradient computation for the synthesis approach is detailed 
in Section 4.4 and that for the integrated synthesis approach 
is given in Section 4.5. 

The fully constminod design formulation is not amenahle 
to such explicit gradient computation procedures as used for the 
other tv/o formulations hccausc it is basically an itorotivo 
design pi’ocodurc. furthermore, the numher of design variables 
in this formulation is considerably small, and any saving of 
computer time obtained in the computation of gradients would 
not bo significant. The gra-dient is, therefore, obtained for 
this formulation by the forvrard difference scheme, ' 

4.5 METHODS OE KEAIALISIS: 

The problem of reanalysis may be stated as follows: 

Given a design of a structure, its stiffness ina.trix, and its 
displacements and stresses under a set of loads, wlmt is the 
Dest method of obtaining the deformations and stresses of a 


modified design? 



49 


The geneml methods of solution adopted for this problem 
are classified into a) Taylor* Series Expansion or the First 
Order Variation Methods, b) Iteration, c) Power Series Expan- 
sion, and d) Reduced Easis Methods. 

The Taylor series expansion methods are more v/ide spread 
than the others. Consider the set of equilibrium equations, 

K u = p (4.3) 

of the starting design. Y^/hen the design is modified, the 
stiffness matrix changes from K to K and the deformations 

change from u to u + ^u. The equilibrium equations of the 
modified design are given by 

[k +Ak] ^u +AU} = 5 (4.4) 

or K u + AK u + K du + AKau = p (4.5) 

By virtue of Equation 4.3 mid ignoring tlio second order term 
AKAu, and rearranging, we obtain, 

Au = - K“^AK u ( 4 . 6 ) 

.-.1 

The analysis of the starting design would have K and u in 
the memory of the computer, or, if the . solution of Equation 4 ., 3 
is obtained using tho Cholesby decomposition, the decomposed 
upper or lower triangular matrix and u would bo available. 
Computation of AK followed by matrix multiplications would- give 
Au and thereby u + Au. This procedure is applicable to small 
modifications in the design in which the effect of ignoring the 
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tem AKAu vrould noi; Idq significant, for largo changes in the 
design, when the Glements in hK would be compamblc to those 
in matrix K, the following equation is used: 

Au = _ [K +AK]~''£iK u (4.7) 

The matrix [^K is computed from IC"'^ using the Sherman 

and Morrison formula j63|. Seveia .1 variations of this basic 
procedure have been used by Kirsch and Rubinstein j 64] , Argyris 
ot al j 65,66,671 and Kavlie and Powell jOsj. 

Kirsch and Rubinstein | 69 | liave considered the problem 
of ronn£:(,lysis by iteration. The iteration scheme is basically 
given by, 

= u + K"'' ZiK u^^^ (4.8) 

in which the u with a superscript denotes the displacement for 
the modified design and the superscript denotes the iteration 
count, for large changes in the design, the convergence of 
the iterative procedure was found to be poor, or oven not . 
possible. 

Romstad et al j 70 l have used a power series expansion 
procedure to compute static and dynamic structural response of 
modified structures. This, the authors state, was necessitated 
because the first order Taylor series expansion was found to be 
inadequa,te for the optimum design v/ith stability constraints. 
Kirsch and Rubinstein's iterative scheme j 69 l given by 
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Eqiiation 4,8 can also bo interpreted as a kind of pov/er series 
expansion. 

The rcdncci basis method was proposed by fox and Mirra, (7lj 
in which the anal37'sts of a modified design is expressed as a 
linear combination of the exact analyses performod earlier for 
different designs, and the problem reduces to tho solution of 
a smaller system of equations. The method lacks tho mathemati- 
cal rigour, and the accuracy of the solution obtained is not 
predictable. Recently, Noor and lowdor j 72 j ha.vo proposed a 
rational choicu of the reduced basis vectors. Ihc method would 
bo competitive only for problems in which the number of degrees 
of freedom is very largo as compared with the nuir.bor of design 
variables. 

Roor and lowder,j72j have also considorod the accuracy 
and convergence characteristics of a rcanalysis procedure com- 
bining a single term laylor series expansion follov/od by itera- 
tion, and have obtained accurate results with vorjr few iterations 
even for significant changes in the design variables. 

In problems of reanalysis, it is very useful to consider 

tho concept of 'Sensitivity Analysis Vectors'. Essentially 

tho first order sensitivity analysis vector for a design 

variable is tho rate of change in tho analysis vector (c.g. 

displacements) with respect to a change in that design variable. 

l.i. T. '■'7^'FUR 

This is also the result of expanding the analysi^pp^-^o^ 

lOT 

Aoc. Mo. aZ t.*v * « 
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laylor scries at a design point with respect to the design 
varia:blGs. Per example, 


u (x +Ax) = u(x) + E |:r ^x. + error (4.9) 

i=:-l ^ 

011 1 

0 ^ 1 ^ is defined as the sensitivity ana, lysis vector for dis- 
placements at point X with respect to the i-th design variable. 
This concept has been used for static roanalysis by Noor and 
Lowder j 72 j and for gradient computation in the present work. 

A similar concept is used for eigenvalue analysis by fox and 
Kapoor j 73i and for flutter ara-lysis by Eao j 74j . 


Anoth>ir procedure v^hich is gaining popularity of late is ■ 
the method of successive approximations. This procedure tries i 
to reduce the total number of analyses performed during an ; 

optimization procedure. The objective fu.nction and the cons- 
tiaints of the optimization problem are expressed in terms of ; 
linear or quadratic Taylor series expansions. The minimization j 
is performed for the approximated function and constraints. 
Successive approximations and minimizations lead the solution 
to the optimum of the actual problem.. The advantage of these ; 

procedures is that very few number of analyses is required to ; 

make the approximations, and the optimization procedure itself I 
does not require any additional analysis. Thus the number of 
analyses performed during an optimization procedure is considera.blj 
reduced. The linearisation concept has been used by ' i 
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Romstad and Fang j 1oi and iDy Sclamit and Farshi | 75j . 'Ihe 
qmdratic expansion techniq-ae has been studied by Lund {76j. 

This approach seems to be very attractive when combined ?/ith 
the Sequential Unconstrained Minimization Technique. 

The procedures used for the computation of gradients 
in the present work are given in the next two sections. In 
the synthesis approach the Taylor series expansion is used to 
obtain the sensitivity analysis vectors for displacements and 
then the finite difference scheme is used to obtain the rates 
of change of the constraints. In the integrated synthesis 
approach, closed form differentiation of the expressions is 
used where available. Finite differences are used where closed 
form differentiation is not possible. In both the cases, the 
unnecessary computations are avoided. 

4.4 GRALIEB'T COl/lTUTATION FOR THE SYNTIIESIS APPROACH: 

Let X be the current design vector, where the gradient 
is required. Let K and u be the stiffness matrix and displace- 
ment vector respectively at this point. A small change f'.x. 

J 

in the d-th design variable causes a change of klvK. in the 

cl 

stiffness matrix and the change in the displacement vector 
is given by, 

Au. = - K“^ Ak. u 

3 3 


(4-. 10) 
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In this formulation, the design variables excluding the 

pretensions in the cables are all properties defining the 

cross section of any one member. The ctuantity fiK. can there- 

J 

fore be computed, for these variables, from the stiffness matrix 
of the element i, to which the variable 3 is related. In the 
linear elastic analysis of the structure, the stiffness matrix 
of the element is linearly related to its area of cross section 
A and moment of inertia I. 


S. 
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(4.12) 


eA. 

^ 6x . 


ai. 


(4. 15) 


Thus AS.., the change in the stiffness matrix of element i 

due to a change Ax. in the variable 3 can be obtained from 

J 

the expressions for the stiffness matrix by simply substituting 

5A. 61. 

9 "~Ax. and fon A^ and respectively. In the non- 

^ j J j J 

linear analysis, the expression for the derivative of the 
stiffness matrix is slightly more involved due to the presence 
of the axial load parameter 0 . Ihe relevant expressions 
requix’Gd in this case are given in Appendix A. 
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Since only ihe stiffness of member i is affected by 
modifying the variable 3 , . alone contributes to the 

It] 

matrix . The vector Au^ is computed using Equ^ition 4 . 10 
and the displacements for the perturbed design are obtained 
as u. = u + du. . 

, tJ D 

If the design variable is the pretension force in the 
cables, only the load vector changes from p to p + v^rhen this 
variable is changed. Ihe displacement corresponding to this 
change is directly obtained. 




u + du. 
J 


u + Ap 


(4. 14) 


For a change Ax. in the design variable x., the 

J J 


constraints at the perturbed design point are evaluated as 
, k = 1, FOG. The- rate of change of the constraints 

is, therefore, 


5x - 
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^ki “ % 
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(4.15) 


The cost of member i is 


C. = c. A. (4.16) 

where c^^ is the cost per unit volume of member i, and Ij^, 
its length. The rate of change of the cost of member i with 
respect to the change in the variable j is 

60 . 61 . 

1 _ ^ T U 

6x . °i "“i 6x. 

0 2 


(4.17) 
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Again since only memler i is related to the variable 3 , this 
expression is also the rate of change of the objective function 
with respect to the cha.nge in the variable j, 


5f(x) _ 

9x. ax. 

J J 


(4.18) 


Substituting the expressions 4.15 a.,nd 4,18 in Equation 4.2 
and noting that no equality constraints are present in this 
formulation, v/e obtain, 




r 


NOG 

S 

k=1 


%c.i 

Aij T'il 


(4.19) 


which is the j-th element of the gra.dient vector of the penalty 
function. 


4.5 G-RADIEIIT GOMUTATION EOR IHE INlE&iUlIBi) SYNIIISSIS APPSDACH: 

The design vector x in this forraulation consists of three 
types of variables: nodal displacements, geometrical proper- 
ties defining the cross sections of members and the pretension 
forces in the cables. These are treated separately for the 
purpose of development of the gradient expressions. 

let X- be the nodal displacement at joint JT. A change 

J 

in this displacement causes a change in the equilibrium state 
of the related joints. The equilibrium constraint is given 
by Equation 2.12 as, 

1 - 

\ C ^kr ^r " ^k ] ’ k = 1 , . . . ,M 
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The memher end forces and consequen'tly the stress 
resultants change only for member i. Therefore, the stress 
constraints at the perturbed design g, . need only to be 
calculated for member i, and the rate of change of these 
constraints is then given by Equation 4*21. Eor the cons- 
tiaints related to all the other members, the ra,te of change 
due to a change in x. is zero. The rate of change of the 

J 

objective function in this case is given by Equation 4 . 18 . 

Einally, let x. be the pretension in the cable i. A 

J 

change x . in the pretension would cause a change Ap in the 

cl 

load vector and the rate of change of the equilibrixan equations 
would then be 


^0 _ .__J 

8 x. IMAXt ax. 

D ^ D 


M (4,24) 


This change also affects the two stress constraints of 
the cable. The rate of change of the two constraints can be 
v/ritton in closed form as, 




X 0.6 < r ^ 


k e hCC 


(4.25) 


1 % 

8x. 


A. X 0.8 cr^ 


and 


k £' NCC 


( 4 . 26 ) 



59 


Eq_'uation 4.25 relates to the constiaint or the protensior 
force and Eq_uation 4.26 to the constra.int on the final force 
in the cahle. The other stress constraints and the objective 
function are unaffected hy a change in this variable. 

'^he gradient of the penalty function with respect to the 
variables is obtained by substituting tho appropriate rates 
of change as detailed above in the ox|jrcssion for the gradient 


giVvjn by Equation 4,2. 
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ILLU STRilTI YE EXMAPLE S 


I'he foCTiulations discussed earlier are applied to the 
mirini-uir cost desigr* of three different cahle stayed structures. 
Furthermore, the synthesis approach is used to study the effect 
of variation of some of the predetermined paramo tors on the 
optiriiuiri design. The res^ilts of these studies arc presented in 
this chapter. All the results given herein were ohtained on 
the IHY 370/155 system. The execution times indicated for 
various prohloms aro the C.P.U. times on this s^rstom. 

5.1 EmiPIE 1 _ I'lEE CAHTIIEYER: 

'the first structure considered is a cantilever box beam 
supported at the tip by a high tensile ca,blc, as shown in 
Fig. 5.1. the relevant details of the structure arc as follow: 

Member 1 is a box beam 15 cm wide and 30 cm deep, made 
up of mild steel of relative cost 1000 units/m^, and member 2 
is a high tensile cable, under initial tension, of relative 
cost 5000 units/m^. The allowable stresses in the members are: 


PjOx boam: Maximum stress in combined bending 

2 

and axial force = 1250 Kg/mn 

Maximum average shear stress in 

= 945 hg/cm^ 


the web 



61 


p 

Cal^le: Yield stress of the cable, CT = 15000 leg/ cm 

? 

MaxtIn^¥n pretension stress = 0.6 CT^ = 9000 kg/cm 

yc ^ 

Maximim stress under working 

load = 0.8 (T = 12000 kg/cm^ 

y c 

Iho vortical deflection at joint 3 is constrained so as not to 
exceed. + 1.5 cm. 

The synthesis approach for this example has 4 design 
variables, 4 stress constraints, 1 deflection constraint and 
4- side constraints. Ihc intogra,ted synthesis approach leads 
to a problem Vi/ith 7 design variables, 3 oq^uilibrium constraints, ■ 

4 stress constraints, and 7 side constraints. Iho fully cons- 
trained design formulation, which loads to an unconstrained 
minimization problem, has only one design variable for this 
example , 

Ihvo typos of loads are considered on this structure 
leading to two loroblcms. In Problem 1, the structure is loaded 
by a single downward load at joint 3 of niagnitudo 40 Tonnes, 
and in Problem 2, the beam is loaded b 3 r a, downv;a,rd uniformly 
distributed load of ST/m. For each of the problems, two cases 
are studied; case a in which the minimum thicknesses of the 
flange a.nd wob of the beam are restricted to 6 mra and case b 
v/horc no such restriction on the thicknesses is imposed. Both line- 
ar and nonlinear analyses arc considered for obtaining the 
optimum design via the synthesis and the integntod synthesis 
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approaches. However, only the linear bcho. viour is considered 
for the fully constrained design formulation, for reasons 
to be explained later. 

Ihe results of the throe formulations for Problem 1 
are presented in Ta.blcs 5.1, 5.2 and 5.3, and those for 
Problem 2, in Tables 5.4, 5.5 and 5.6. It may bo observed 
that with the exception of Problem 1 case b, all the throe 
formulations lead to the same design in all the problems (within 
a certain nimcrical ccccuracy) and also that the final design 
is a fully constrained design. It may also bo observed that 
the optimum design in tho ease of Prohlcm 1 ease h is such 
that the area of the box beam is nearly constant although the 
thicknessGs of tho v/ob and flange arc quite , different in 
all the cases. This is discussed in greater detail in 
Section 6,2. 

5.2 EX/iiiPIE 2 - SmiETRIG CABLE STATED &IRDER: 

The second structure is a part of the symmetric cable 
stayed girder as shown in Pigure 5.2. Only the portion of the 
girder to the left or right of the suspended span is considered 
for optimization. This portion is shown in figure 5.3. 

The girder is a box beam 200 cm wide and 300 cm deep, 
and is supported by 3 high tensile cables under initial tension. 
The cables meet at the top of the tower, and the tower is also 
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assumed as a box section 200 cm wide and 300 cm deep. The allow- 
able stresses in the beams and cables are the same as those used 
for the previous example. 

The synthesis approach to the optimviin design of this 
structuro leads to a nonlinear programming problem with 16 design 
variables, l6 stress constraints and 16 side constraints. The 
integmted synthesis approach has 31 design variables, 15 equili- 
brium constraints, l6 stress constraints and 31 side constraints. 
The fully constrained design formulation has only three design 
variables. 

Once again, two types of loads, (i) only joints loaded 
and (ii) members and joints loaded are studied for this structure. 
Problem 3 considers the total loads on the structure to be lumped 
at the nodes of the structure; thus joints 2,5,5 and 6 are sub- 
jected. to downward loads of 3851', 6oT, 340T and 330T respectively. 
Problem 4 considers the loads on the girder to be uniformly 
distributed at loT/m. Besides, joints 3 acad 6 have downward 
loads of 6oT and 300T respectively. 

The integrated synthesis approach failed to give meaningful 
results for this problem because of premature termination of the 
program. The experiences v\fith this formulation are discussed in 
the next chapter. Therefore, the results are presented only for 
the synthesis approach and for the fully constrained design 
fo rmula tion. 
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The results, using the synthesis approach, of the linear 
and nonlinear formulations for Problera 3, ?vdth costs of the 
materials as 1 and 3 units/m^ for beaxi and cable, are presented 
in Table 5.7. The results show that the minimum cost design 
is a fully constrained design. To ascertain this fact, a few 
variations are studied. If the design is a fully constradned 
design, the variation in the costs of the materials should not 
affect the tota.1 quantity of the materials used. Indeed it is 
found to be so, when the relative cost of the cable is changed 
from 3 to 5 units/m^. This result is x>nesented in Table 5.8. 

The above cases consider the lov/er bound on the thickness 
of the flanges and the webs to be 16 rora. This bound is reduced 
to 10 mm and the problera is solved again for the two cost ratios 
mentioned above. These results are also presented in Table 5.8. 
All these results confirm that the optimum design is a fully 
constrained design. 

It may be observed that for this problera, the nonlinear 
behaviour has very little effect on the optimum total cost. 

The optimum total cost of materials is roducod slightly if the 
nonlinear behaviour is considered. 

The results, using the synthesis approach, under linear 
and nonlinear behaviour, of Problem 4 are presented in Table 5.9. 
The rfliniraum thickness limitation of l6 mm and tho unit costs of 
1 and 3 units por m^ for beam and cable respectively are 
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considered, The reduction of the cost of the structure with 
successive values of the penalty pammetor 'r* a.re shown in 
Figure 5.4. The effect of nonlinearity in this prohlem is 
once again ohservod to be very small, Hovrovor, the consideration 
of nonlinear effects results in an increase in the optimum total 
cost of the materials. 

The results of Problems 3 and 4 using tho fully cons- 
trained design approach are presented in Table 5. 10. The minimum 

thickness limitation is 16 mm and tho unit costs arc 1 and 3 units 

■5 

per m-^ for tho beam and tho cable respectively, lor Problem 4> 
the reduction in the cost of the structure v\rith respect to the 
cumulative number of linear searches performed for unconstrained 
minimization is shown in Figure 5.5. Tho variation of tho design 
variables during the search procedure is shown in Figure 5.6, 

It is seen from Figure 5.5 that the fully constrained design 
at tho starting point is only about 1 percent costlier than the 
optimum design. This indicates that the region over which tho 
optimisation is to bo performed is rather 'shallow', 

5.3 EK/dvlPIE 3 _ lOHSIMIfflTRIG CAELE STATED GIEIER: 

The third structure considered is shown in Figure 5.7. 

It is a slight variant of the one studied by Tang | 25j and 
Morris 127]. The structure consists of 4 beam elements, 1 tower 
and 2 cables under initial tension. The beams and the tower 
are box girders 250 cm wide and 350 cm deep. The structure is 
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loaded 03^ a uniformly distributed load of 25iym and tlie top 
of the tower is subjected to a downwcird load of 85'!'. Ihe lower 
bound on the thickness of flange and web of the beam elements 
is 12 mm. The allowable stresses are the same as the ones 
used in the earlier problems. This problem is henceforth 
referred to as Problem 5. 

The synthesis approach as applied to this problem has 
14- design variables, 14 stress constra.ints and 14 side cons- 
traints. The results of optimization for the linear and non- 
linear behaviour of the structure are presented in Table 5.11. 
It may bo observed again that the effect of considering the 
nonlinear behaviour on the optimum design is very small, and 
that the optimum design ha,s a slightly increased cost. The 
reduction in the cost with respect to successive values of the 
penalty parameter is shown in Figure 5.8. 

The fully constrained design formulation as applied to 
this problem has only two design variables. The results of the 
problem using this formulation are presented in Table 5.12. The 
reduction in the cost of the structure v/ith rospoct to the 
cumulative number of linear searches performed is shown in 
Figure 5.9j and the variation of the design variables during 
the search procedure is shown in Figure 5.10- It is seen 
from Figure 5.9 that the fully constrained design at the start- 
ing point is 4.5 percent costlier than the optimum design. It 
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may be observed that tbe time taken for arriving at the optimum 
is only about a tenth of the time taken by the synthesis 
approach. 

In order to study the effect of the depth of the section 
on the extent of nonlinearity that would be present in the 
optimum design, a pa,ranietric study of the variation of the 
optimum design with respect to the depth of the girder is made. 

The girder is assumed to be a box 600 cm wide, and the depthx is 
varied from 150 cm to 350 cm at intervals of 50 cm. Ihis problem 
is referred to as Problem 6. The results of these studies for 
the linear behaviour of the structure are presented in Table 5.13 
and for the nonlinear behaviour, in Table 5.14-. fhe variation 
of the optimum cost with respect to the depth of the girder is 
shov/n in figure 5.11. It may be observed that the effect of 
nonlinearity on the optimum design is only margins,! in all cases ^ 
C'-nd t'hat the pretensions in tbe cables undergo relatively mop^' ■ 
significant changes for smaller depths. 

The structure is optimised for a depth of girder of 

150 cm and v/idth of 600 cm, v/ith the yield stress '.rg of the 

yC 

O O 

cable reduced to 10000 kg/cm'^ from 15000 kg/cm'^ considored 
earlier. The results for both the linear and nonlinear analyses 
are presented in Table 5.15. It is scon that the total quantity 
of mild steel remains pra,ctically unaltered and that the quantity 
of high tensile cable increases almost in proportion to the 
reduction in the allowable stresses. 
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box BEAM 


X 1 ; 




_ 500 cm 


LOADING: . 

PROBLEM 1 ■, DOWNWARD LOAD 40 TONNES AT JOINT 3 

PROBLEM 2: DOWNWARD U.D-L. B TONNES /m ON MEMBER 
FIO. 5-1 CANTILEVER BOX BEAM SUff^ORTED AT THE TIL 
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Table 5.' 

1; Results of Problem 1 
Approach. 

by the Synthesis 

Member Qty. 

Linear 

Nonlinear 

Variable . 

Constraint 

Variable 

Stress 

Constraint 


Case a minimum thickness 

= 0.6 cm 


1 t , cm 

w’ 

0.6000 0.0001 

0.6000 

0.0000 

■ t^, cm 

0.6000 0.9863 

0.6002 

0.9863 

2 A , cm^ 

4.8573 0.0001 

4.7971 

0.0000 

P, Tonnes 

43.714 0.0409 

43.179 

0.0288 

Tef . Constr . 

0.1663 

0.0804 

Cosb 

44.1733 

43.9639 

bo. of functions 

230 

253 


No, of Gradients 

61 

50 



Case b minimum thickness 

= 0.0 cm. 


1 t , cm 

0.1881 0.0000 

0. 5903 

0.0002 

t^, cm 

0 . 7980 0.9909 

0.0234 

0.9984 

2 

2 A , cm 

6.2117 0.0000 

6.0397 

0.0002 

P , Tonnes 

55.905 0.2430 

54. 355 

0.2206 

Def, Constr. 

0.9572 

0.9171 

Go st 

39.5748 

39.4142 

No. of functions 

257 

279 


No. of Gradients 

50 

63 


Note: The non-dimensional stress constraints in tbe tables 

are in tMe order given in TablG^2.1 for the type of 
member. 
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Table 5. 3 1 Results of Problem 1 by the Pully 
Constiained Pes.lgn Approach. 


Member 

Qty, 

Case a 

Case b 

1 

t , cm 

w’ 

^f 

0.6000 

0. 6001 

0.0024 

1.1795 

2. 

p 

A, cm 

P , Tonnes 

4.8569 

45.712 

6.1612 

55. 450 

Starting 

Design 

Pg, Tonnes 

40.0 

40.0 

Cost 


44.1729 

39.5479 

No. of Designs 

119 

233 


Note: The number of designs performed is used as an 
indication of time rather than the number of 
function and gradient evaluations. 



Table 5.4: Results of Problem 2 by the Synthesis 
Approach. 


Member 

Qty. 

linear 

Nonlinear 

Vari- 

able 

Const. 

Vari- 

able 

Const. 




Case a 



1 

t , cm 

w’ 

0.6011 

0.0000 

0.6002 

0.0002 


t^jcm 

3.0073 

0.3125 

3.0087 

0.3114 

2 

2 

A , cm 

3.8460 

OiOOOl 

3,8722 

0.0002 


P, Tonnes 

34.611 

0.4923 

34.843 

0.4957 

Pef. Const. 


0.1874 


0.1636 

Cost 


76. 

7392 

76 

.8252 

No . o f 

Functions 

269 


323 

No . of 

Gradients 


70 


80 




Case b 



1 

t , cm 

0.4134 

0.0001 

0.4133 

0.0000 


t^jcm 

3.1572 

0.0003 

3.1567 

0.0000 

2 

p 

A , cm 

3.8327 

0.0001 

3.8591 

0.0000 


P, Tonnes 

34.492 

0.4906 

34.731 

0.4941 


Def. Const. 0.2001 0.1761 


73.3098 

334 

77 


75.3929 

476 

97 


Go st 

No. of Functions 
No. of Gradients 
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Talole 5,5- Results of Problem 2 by the integrated 
synthesis approach. 


Member Qty.^ 

LINEAR 

NONLINEAR 

VAR 

CONS 

VAR 

CONS 

, cm 

-0.0312 

Case a 
-0.0009 



, cm 

1.3522 

-0.0004 



U3 

0.0084 

0.0013 



t , cm 

1 

t^, cm 

0. 6004 

3.0086 

0.0002 

0.3117 

PlffildATURE 

TEIBI NATION 

0 

A ,cm 

3.8472 

0.0005 



^ P , To nne s 

34-608 

0.4925 



Cost 

76.7424 



No. of functions 

436 



No. of Gradients 

112 





Case b 



IX..J 9 cm 

-0.0329 

-0.0041 

- 0.0329 

0.0085 

, cm 

1 - 3409 

-0.0224 

1.3558 

0.0090 


0.0083 

0.0133 

0.0084 

-0.0008 

t , cm 
. y-r 

t^, cm 

0.4147 

3.1579 

0.0012 

0.0035 

0.4196 

3.1586 

0.0012 

0.0156 

p 

A, cm 

2 

P, Tonnes 

3.8396 

0.0024 

3.8798 

0.0050 

34.474 

0.4922 

34.744 

0. 4968 

Cost 

73.3825 

73 

. 6865 

No. of Functions 

387 


412 

No. of Gradients 

89 


95 
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Dalole 5.6: Results of Problem 2 by the Pully 
Constrained Design Approach. 


Membe r 

Qty. 

Case a 

Case b 

1 


0,6000 

0.4133 


t^, cm 

3.0096 

3.1564 

2 

P 

A , cm 

3.8421 

3.8337 


P , Tonnes 

34.579 

34.503 

Starting 

Design 

Pg, Tonnes 

40 '.0 

40.0 

Cost 


76.7279 

73.2982 

I'lo. of Designs 

87 

97 
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Table 5.7: Results of Problem 5 by the- Synthosis 
Approa ch 

Cost Ratio 1:3; Minimum Thickness 1.6 cm. 




Initial 

Bosisn 

linear 

iioni inear 

Mem- 

ber 

Qty 

Vari- 

able 

Stress 

Cons- 

traint 

Vari- 

able 

Stress 

Cons- 

traint 

Vari- 

able 

Stress 

Cons- 

traint 

1 

t , cm 

w’ 

2.0 

0.5371 

1 . 60 1 6 

0.3371 

1.6023 

0.3360 


^f 

4.0 

0.9684 

1.6064 

0.9791 

1 . 6061 

0.9790 

2 

t 

w 

2,0 

0. 3284 

1 . 6007 

0.0000 

1 . 6027 

0.0000 


t. 

X 

4.0 

0.9847 

1.6037 

1.0000 

1.6055 

1.0000 

3 

t 

w 

2.0 

0.0561 

1 .6005 

0.0001 

1.6031 

0.0007 


t^ 

4.0 

0.9697 

3.9757 

0.9783 

3.9809 

0.9768 

4 


2.0 

0. 4720 

1,6119 

0.0022 

1.6525 

0.0231 


^f 

4.0 

0. 7090 

1.8405 

0. 6389 

1.8716 

0.6477 

5 

t 

w 

2.0 

0. 5102 

1.6007 

0.1384 

1.6016 

0.1874 



4.0 

0.9927 

1.6032 

0.9909 

1.6062 

0.9941 

6 

2 

ii. , cm 

150.0 

0.4444 

77.180 

0.0038 

74 . 63 s 

0.0061 


Tonnes 

750.0 

0 . 5268 

691.98 

0. 1522 

667 . 63 

0.1232 

7 

A 

280.0 

0. 4444 

131.98 

0.0139 

132.40 

0.0017 


P 

1400.0 

0.5187 

1171.31 

0.0002 

1 1 89 . 62 

0.0014 

8 

A 

125.0 

0.4444 

75.608 

0.0008 

73.266 

0.0005 


P 

625.0 

0.4870 

679.97 

0.1006 

659.08 

0.0721 

Yolume of 
Beam (m3) 

49. 

000 

33 

.980 

34 

.040 

Volume of 
Gable (m3) 

3. 

836 

1 

.984 

1 

.952 

Cost 


60.5084 

39 

.9312 

39 

.89 79 

No . 0 f 
Functions 




450 


499 

No . 0 f 
Gradients 




126 


136 

Time 

, Secs. 




254 


336 
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Table 5.8: Results of Problem 5 by the Synthesis Approach. 

(Influence of Cost ratio and minimum thickness). 


Cost Ratio 1:5 Cost Ratio 1:5 CSost Ratio 1:5 
Min.Thk. 1.0 cm Mi n.Thk.1.0 cm Mi n . T hk.1.6 cm 


Mem- 

ber 

Qty. 

Vari- 

able 

Stress 

Const. 

Vari- 

able 

Stress 

Const. 

Vari- 

able 

Stress 

Const. 

1 

t , cm 

w’ 

1.0016 

0.0792 

1.0007 

0.0011 

1.6079 

0.3183 


^f 

1.0032 

0.9859 

1.0051 

0.9715 

1.6070 

0.9776 

2 


1.0468 

O.OOOQ 

1.0394 

0.0001 

1.6077 0.0001 


^f 

1.7383 

1 .0000 

2. 1276 

0 .9999 

1 . 6680 

1.0000 

3 


1.0034 

0.0002 

1.0007 

0.0002 

1.6021 

0 . 0013 . 


■^f 

4. 5479 

0.9499 

4 . 5368 

0.9562 

3.9853 

0.9768 

4 

"t 

w 

1.0074 

0.0108 

1.0045 

0.0048 

I.S 311 

0.0052 


h 

2.1652 

0.4222 

2.1507 

0.4205 

1,8385 

0.6431 

5 

t 

w 

1.0213 

0.0008 

1.1189 

0.0008 

1 . 60 60 

0.2101 


% 

1.3918 

0.9999 

1.2349 

1.0CC0 

1. 6061 

0.9956 

6 

2 

A, cm 

84.505 

0.0005 

78.183 

0.0010 

74.421 

0.0019 

P 

, Tonnes 

760.18 

0.2468 

702.95 

0.1690 

668.52 

0.1176 

7 

A 

153.88 

0.0001 

147.56 

0.0007 

135.44 

0.0009 


P 

1384.73 

0.1516 

1325.36 

0.0978 

1217.91 

0.0238 

8 

A 

90.447 

0.0012 

78.542 

0.0007 

72.359 

0.0016 


P 

813.01 

0.2203 

706.35 

0.1260 

654.65 0.0713 

Volume 
Beam ( 

of 

m3) 

28 

.726 

29 

.427 

34 

. 190 

Volume 

Cable 

of 

(m3) 

2 

.300 

2 

.120 

1 

.969 

Cost 


35 

. 6256 

40 

.0267 , 

44 

.0344 


No. of functions 472 614 426 

No. of Gradients 120 I 68 ; 114 
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Table 5.9: Results of Problem 4 by the Synthesis Approach. 
Cost Ratio 1J'5; Minimum Thickness 1.6 cm. 




Initial Design 

Lines r 

Nonlinear 

Mem- 

ber 

Q'ty. 

Vari- Stress 

able Const. 

Vari- 
able • 

Stress 

Const. 

Vari- 

able 

Stress 

Const. 

1 

t - cm 

w’ 

2.0 0.3566 

1.6006 

0.0004 

1 . 6005 

0.0006 



4.0 0.8247 

1 . 60 2 6 

0.8034 

1.6116 

0.8031 

2 


2.0 0.0572 

1 . 6002 

0.0002 

1 . 6007 

0.0002 


"f 

6.0 0.7501 

3.9510 

0. 7060 

3.9545 

0.7062 

3 

t 

V.' 

2.0 0.0814 

1.6005 

0.0002 

1. 6005 

0.0001 


■^f 

7.0 0.6989 

4.7131 

0. 6483 

4.7133 

0.6483 

4 


2.0 0.4720 

1 .6037 

0.0025 

1. 6071 

0.0065 


tf 

4.0 0.6825 

1 .8442 

0. 6042 

1.8534 

0.6052 

5 


2.0 0.4663 

1 . 6009 

0.0004 

1 . 6009 

0.0014 


h 

4.0 0.9842 

1 . 60 1 1 

0.9803 

1. 600 7 

0.9804 

6 

2 

A ,cm 

150.0 0.444^ 

125.76 

0.0001 

127.49 

0.0001 

]? 

, Tonnes 

750.0 0.5515 

1131.80 

0. 4697 

1147.29 

0.4766 

7 

A 

280.0 0.4444 

206 . 55 

0.0000 

210.95 

0.0001 


P 

1400.0 0.5636 

1858.89 

0. 3878 

1898.42 

0. 4006 

8 

A 

125.0 0.4444 

113 . 41 

0.0000 

114. 68 

0.0002 


P 

625.0 0.5382 

1020. 62 

0. 4500 

1031.97 

0.4566 

Volume of 
Beam (in3) 

60 . 440 

40 

. 49 1 

40.513 

Volujae of 
Cable (m3) 

3.836 

3 

.093 


3. 144 

Cost 


71.9484 

49 

.7710 

49.9453 . 

ITo . 0 f 
functions 



524 


570 

lo. of 
Gradients 



143 


159 

Time 

, Secs. 



289 


414 
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Tatle 5.10: Results of Problems 3 and 4 "by the Pully 
Constrained Design Approach. 


Member 

Qty. 

Problem 3 

Problem 4 

1 


1.6000 

1 . 6000 


^f 

1,6000 

1.6435 

2 


1 . 6000 

1.6000 


^f 

1.6417 

3.9839 

3 

t 

w 

1 . 6000 

1 . 6000 


h 

4.0012 

4.7212 

4 

t 

w 

1.6000 

1.6000 


^f 

1. 8402 

1.8398 

5 


1 . 6000 

1 . 6000 


^f 

1 . 6000 

1 . 6000 

6 

A, cm^ 

70.731 

123.28 


P, Tonnes 

636,58 

1109.52 

7 

A 

153.55 

208.93 


P 

1381.98 

1880.39 

8 

A 

73.405 

110.46 


P 

6 6o . 6 5 

994.12 


Pg, Tonnes 

1 

'-^8 

750.00 

1000.0 

Starting 

1400,00 

2000.0 

Design 

625.00 

1000.0 

Cost 


40.2427 

49.8249 

ho, of Designs 

108 

118 

'Tiine 5 Secs. 


32 

36.3 






FULIY CONSTRAINED DESIGN 
AT STARTING POINT 


NUMBER OF LINEAR SEARCHES 


FIG.5'5 PROGRESS OF OPTIMIZATIGN-FULLY CONSTRAINED 
' DE’SIGN APPROACH/FOR PROBLEM 4 






MEMBER 7 


MEMBER 6 


MEMBER 8 


NUMBER OF LINEAR SEARCHES 


3.S-6 PATft OF THE DESIGN VARIABtES-FyLLY: eONSTRAlNED 
OFSifiN APPROACH FOR PROBLEM A 
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'•’’aljlc 5, 

11: Results of Problem 5 by the Synthesis 
Approach. 


-uiii- 4ty 

* t'; V 

Vari~ Stress 
able Const. 

__...Pinear 

NonlineaT 

vari- 

able 

Stress 

Const. 

Vari- 

able 

Stress 

Const. 

1 

«L 

2.0 

0.3160 

1.3266 

0.0002 

1.3270 

0.0000 

y « 

1 

10.0 

0.3301 

6.4713 

0.0002 

6.4595 

0.0000 

f'.’ "li 

..i. 

2.0 

0.1596 

1 , 2009 

0.0000 

1.2003 

0.0000 

1^' ^ p 

15.0 

0.4218 

12. 4190 

0.0559 

12.4487 

0.0554 

5 t,. 

# 1. 

.f:.. 

2.0 

0.3126 

1.2021 

0.0001 

1.2005 

0.0000 


15.0 

0.6356 

12. 4160 

0.3039 

12.4422 

0.3019 

^■f t 

1 # 

2.0 

0.5968 

1.2006 

0.0003 

1.2002 

0.0001 


4'» 0 

0.6617 

1.5278 

0.4320 

1.5269 

0.4318 

'i t 

V.’ 

4.0 

0.0763 

1.2689 

C .0006 

2.9974 

0.0000 

'®^:r 

■f"'! 

6.0 

0.9702 

6.1782 

0.9999 

3.7635 

0.9998 

w 4 , cnr* 

1000.0 

0.5000 

521.40 ■ 

0.0001 

521.71 

0.0000 

I' 

4500.0 

0.6295 

4692.14 

0.2997 

4695 . 26 

0.2999 

f ii, 

1000.0 

0.3444 

541.47 

0.0003 

541.44 

0.0000 

T 

5900.0 

0.5510 

4871.81 

0. 1982 

4872.86 

0.1978 

’•nlij:..c' of 








142, 

.620 

102.793 

102.884 

' ' * ,1. ' O f 







'.'ai'iio* (::!'■') 

10. 

.304 

P 

j. 450 

C 

3.451 


173, 

.5318 

119.1422 

119.2408 

• ■ 1 ' , of i4iiictiono 


523 

525 

• H*. uf iiiadii 

unts 


134 

147 

* i #.<* i * ii » ^ IK 



226 

289 
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LINEAR 


CUMULATIVE NUMBER OF LINEAR 
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ANALYSIS 
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lable 5.1S: Results of Problem 5 by the Pully 
Constrained Design Approach. 


o:r. b.;j,r 

Qty. 

Value 

1 

cm 

1.3281 



6 . 4430 

•i 

r 

t 

w 

1 . 2000 



12.4978 

3 

t 

w 

1 . 2000 



12.3236 


t 

w 

1 , 2000 


't.p 

1.5311 

3 

"^w 

1 . 2000 



6.3405 

6 

A,cm^ 

522.22 


P, Tonnes 

4699.95 

7 

A 

544.75 


P 

4902.79 

Dl'arting 

'^6 

4500.00 

Jenign 

“^7 

5900.00 

>^0 1 > t 


119.2891 

^•o. of Designs 
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Diiae, Secs. 


21.22 
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POThl'^iTi 6 _ Effect of Variation of Depth of 
iiricr - Linear Analysis. 


'.i.V 

tj ^ 150 

D =: 200 

D = 250 

D = 300 

D = 350 

*1* 1 

3. 0835 

2.3143 

1.8585 

1.5461 

1.3273 


6.7112 

4.9654 

3.9044 

3.1969 

2.6922 


PJ>4P>6 

1.9 886 

1.5930 

1.3230 

1.2019 

i 

■ f 

H.'>720 

7.3020 

6.3299 

5.6792 

5.1835 

1. 

v: 

1 . A 5 Cl 3 

1.4672 

1.2107 

1.2015 

1-2030 

t. 

'-,.41 5B 

7.3959 

6.4277 

5.7106 

5.1831 

t 

w 

1 . Ls’.'jP 

1.2046 

1-2097 

1.2007 

1.2010 

tf 

1 .blPO 

1.3371 

1.2017 

1.2002 

1.2004 


2.21116 

3.7747 

2.4439 

2.6109 

2.0193 

t* 

3.5046 

2.0514 

2.2956 

2.0132 

2.1365 

? 

rf. j 

372.57 

520.20 

519.13 

■ 519.16 

520.91 

,, « u 1 ; i 4» ; : *f 

4239. 1 

4681.7 

4671.0 

4672.3 

4687.3 



459.29 

496.45 

523.27 

543.56 

1 

4146.9 

4130,8 

4460 - 4 

4708.2 

4888.3 


IP 9. 339 

158.207 

140.467 

129.094 

121.264 

V. 

17'-. 865 

142.829 

124. 681 

112.998 

104.903 

^ 

(:;v^ 

4.825 

5.126 

5.262 

5.365 

5.454 
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Table 5.14: Problem 6 - Effect of Variation of Depth of 
Girder - Nonlinear Ana, lysis. 


Member 

Qty. 

D = 150 

D = 200 

D = 250 

D = 300 

D = 350 

1 

t , cm 
w ’ 

5.0920 

2.3207 

1.9320 

1.5594 

1.3305 


^f 

6. 6668 

4.9407 

3.3849 

3.1914 

2.6846 

2 


2,6550 

1.9908 

1.6730 

1.3317 

1.2053 


^f 

9.0810 

7.3634 

6.3314 

5. 6887 

5.2013 

3 


1.8874 

1.4402 

1.7530 

1.9153 

1.2127 


■^f 

8.2920 

7.2760 

6.3176 

5. 4684 

5.1259 

A 


1 . 5879 

1.2073 

1.3792 

1.2061 

1.2053 : 


^f 

1.8116 

1.3364 

1.2028 

1.2004 

1.2021 

5 


4. 1632 

3.6720 

3.8617 

3.2901 

2.5753 


^f 

3. 3668 

2.3447 

1.7104 

1.7365 

1.9301 i 

6 

1 , cm2 

457,27 

505.99 

521.98 

517.80 

516.51 J 


P , Tonnes 

4114.5 

4552.7 

4694.4 

4660.0 

4647.8 j 

7 

A 

481 .01 

475.71 

492.76 

524.52 

551.08 ; 


3? 

4326.8 

4275.0 

4434. 5 

4720.6 

4653.1 

Cost 


190.332 

158. 692 

141. 502 

129.792 

121.537 1 

VolilEj 

of 






ruami (: 

ffi3) 

175.923 

143.399 

125.704 

113.709 

105.173 

Vo 1'i.Uiio 

0 f 






Gablj 

(r'.3) 

4. 803 

5.097 

5. 266 

5. 361 

5. 455 I 



J in cm. 



350 


bO 250 

DEPTH OF GIRDER, cm 

liON.OF MINIMUM COST WITH RESPECT TO I 
J OF GIRDER- nKOBLEM 6. ■ | 
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Ta-ble 5.15: Problem 6 - Results of Lov^^er Gable 
Stresses 

P = 150 cm and .CJl = 10000 kg/cm^ 

y c 


Member 

Qty- 

linear 

Monlinear 


Analysis 

Analysis 

1 

t , cm 
w’ 

3.0878 

5.0896 



6. 6962 

6*6685 

2 


2.6489 

2.6513 



9.0141 

9.0772 

5 

\ 

1.8017 

1.8718 



8.0130 

8.3410 

4 

\ 

1.5871 

1.5886 



1.8329 

1.8101 

5 


1.2704 

8.7356 


^f 

4.2073 

2.3993 

6 

2 

A , cm 

648. 18 

651.39 


P 5 Tonnes 

5888.5 

3908.2 

7 

A 

687.73 

■ 676.84 


P 

4124.2 

4060. 1 

Cost 


195.458 

196.990 

Volume 

of Beam (m^) 

174. 96 S 

176.561 

Volume 

of Gable (m^) 

6.831 

6.810 



CmPTER SIX 


DISCUSSIONS AND CONCLUSIONS 


IliG results of the example problems proscuted in tho 
previous chapter arc discussed heroin. The outcome of tho 
numerical studies is described and conclusions dravm there- 
from are enumerated. Finally, the suggestions to improve 
upon the techniques used in the present work and to make the 
formulations still more practically oriented are indicated. 

6.1 C0I®ARIS0N OF IHE FOffi/iULAIIONS; 

I'he computational experience shows that tho synthesis 
approach to the optimum design pf the class of cable stayed 
structures studied is quite efficient and reliable. The fully 
constrained design formulation, as proposed in tho present work, 
is certainly more efficient in the sense that, it consumes only 
a fraction of the computer time taken by the synthesis approach 
to arrive at the optimum design. This is due to tho fact that 
only one unconstrained minimization is needed in the case of 
the fully constiainod design formulation. Furthermore, the 
number of design variables in this formulation is relatively 
small. This fact illustrates the point that tho mathematical 
progra.mming methods coupled with an optimality criterion based 
approach would be much faster than the simple use of mathematical 
programming techniques where an optimality criterion can be 
obviously identified. 
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The integrated synthesis approach has more ntmiher of 
design variables and constraints than the synthesis approach. 
This fact itself cancels out the advantages gained in not 
solving a set of simultaneous equations at every design point 
of the search procedure and the use of the much faster gradient; 
computation procedure. The computational time is therefore 
always more than that for the synthesis approach. 

The integrated synthesis approach successfully led to 
the optimal design for problems 1 and 2 hut resulted in a pre- i 
mature termination of the program for one case in Trohlem 2 
and for Problems 3j 4 and 5. This is due to the fact that 
during the first two or three unconstrained minimizations, the ^ 
unbalanced forces and moments in the equilibrium equations, 
which actually form the equality constraint values are brought ^ 
down frcoii'i largo values to almost zero. This is brought about 
without much regard to the objective function, which is the ^ 
cr).‘:;t of the structure. It is found tha.t once the equilibrium | 
c^'nfjtmints , which are nonlinear equality constraints, are 
on tie fled to such a dogreo , tho path of optimization is along | 
tM.,j curved iritorGOction of tho equality constraint surfaces. ; 
Tracing this curved path by straight move vectors loads to ; 
tl’.c pruinature termination of the algorithm. Such a strict ; 
satisfaction of the equilibrium constraints at the early stages 
f th.o optiiiiization procedure is due to the fact that the eon 
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about an order of magnitude less than tho starting values. 

This scheme turned out to be satisfactory for problems 1 and 2 ^ 
only. 

6.2 RESULTS 01 THS IIICSAR AIUiYSIS: 

'the results based on linear analysis show that the 
optimum design for the example problems always results in a 
fully constrained design. This is characterised by the fact 
that , 

a) the webs of the beams are always constrained by the 
maximum thickness limitation or by the shear stress 
limitation; 

b) the flanges of the beams are alvays constrained by the 

iii£!Xiinum extreme fibre stress in combined bending and 
thrust or by the minimum thickness limitation, and ; 

c) the areas of the cables are always constrained by the 
pretension stress or the final stress under working load. : 

The prestressing forces in the cables adjust themselves : 
in such a way that the cost of the resulting design is minimum. ; 
b-iic chaCTcteristic behaviour of this class of cable stayed I 
;,;tr>ict\}rou , in fact, led to the fully constrained design for- | 
mulation used in tho present work. I 

Since the optimum design is a fully constrained design, | 
it is natural to expect that the variations in the relative cos' 
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of mild steel and high tensile cable should not have an : 

influonco on the total quantities of mild steel and cable, 
ihe rcGUlts of Table 5.8 show that this is true. The changes 
in the quantities of mild steel and high tensile cable are only- 
marginal and can bo said to be negligible keeping the large 
variation in the relative costs from 1:3 to 1:5 in mind. 

In case of structures with only nodal loads, the pre- 
tension forces adjust themselves in such a way that at the ! 
optimum design, the beam elements undergo very little bending j 
and aro mainly axially loaded. This is consistent with what ; 
could bo intuitively expected. However, for structures with 
member loads no such conclusion is valid. The comparison of 
the optimum designs for the joint loaded and the member loaded 
cable stayed structures show an absolute lack of correlation, 
and the inomber loaded structure has a much higher cost. It 
is r;uite clear that lumping the loads at the nodes is not 
dosirablo, even in the analysis probloias, and for the purposes | 
of optirnum design, the results of such lumping would lead to ! 
an unsafe design. 

The tower mumbor, at the optimum, is subjected only to | 
on aria] force and the bonding is negligible. Therefore, 
oriLy the area of cross section of the tower is of importance | 

anj. its moment of inertia is only of secondary importance. The: 
optlmu::; dv^sign therefore does not show the typical characterist. 
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of web thickness being at the minimuia and the flange thickness 
being determined by the stress constraint. Any conibina,tion 
of the web and flange thicknesses such that the area of cix)ss 
section is the recj_uired one produces a near-optinimii design. 

In the synthesis approach, convergence may take 'place to any 
point in this plateau. This is the reason for the cross, 
section of the beam in Problem 1 case b., and the cross section 
of the tower in Problems 3 to 6 being different in different 
methods. The area of cross section, however, may be observed 
to be almost the same in all the methods for each of the 
problems . 

Prom the parametric study presented in Table 5.13 it is ; 
observed that the total cost of materials increases as the ! 
depth of the girder decreases. This increase is mainly duo 
to the increase in the qiiantity of mild steel used. The q.uan_ | 
tity of high tensile cable used is seen to decroasc, although | 
only slightly, with the decrease in depth. The pretension 
forces do not show a spocific trend, and their variation is 
rather small for the largo variation in the depth. It is [ 

soon from cxporionco with the fully constrained design formu- ; 
lation that the pretension forces, in the neighbourhood of the 
optiuimi df‘sign, do not significantly affect the cost and that : 
there is a rather shallow region around the optimum. Therefore,; 
a fully constrained design with pretension forces in a small 
range about the optimum would produce a near-optimurn design, : 
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It is for this reason tha-t no specific trend v/ould he observed 
in the values of the pretension forces. 

Comparison of results presented in Tables 5.13 and 5.15 
show that a reduction of the allowable stresses in the cable 

-I 

ty 53-^ percent does not significantly affect the total quantity 
of mild steel used. The quantity of cable used increases 
nearly proportionally, and there is a slight readjustment of 
the pretension forces in the cables. 

6.3 RESULTS Of TilB NOmNEAR ANALYSIS: 

Comparison of the results of the example problems using 
the linear and the nonlinear analysis procedures indicates 
only a marginal difference between tho costs of the two optimm 
deoiipio. However, it is consistently observed that considering : 
the nonlinear effects in a joint loaded structure results in a 
reduction in the minimum cost, v/hile in a structure with member i 
load;:' an increase in the minimum cost is experienced. This is ; 
oxpl;.'incd as follows: ^ 

In a joint loaded structure, the pretension forces at 
the optimum design arc determined such that the structure is 
cs!m,;ntially axially loaded and the beioding effects are small. 

th'j beam-column effect is considered, the flexural stiffnesi 
of the members is reduced, resulting in smaller design moments. I 
The rerjulting design is therefore loss costly. In a structure 
7 ;ii;. fueiaber loads, tho negative moments of the membors govern 
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the design. The negative moments are mainly contrihuted by 
the fixed end moments due to the member loads, and the inclusion 
of beam-column effect increases the fixed end moments. Thus 
the design moments are larger, resulting in a costlier 
structure. 

The influence of the nonlinear analysis of the structure 
on the optimum design is to alter very slightly the thicknesses 
of the 'flanges and the v/ebs of the beam elements of the corres— : 
ponding optimum design using the linear analysis. The nonli- 
nearity effects, it seems, are counter balanced by a small 
variation in the pretension forces in the differont cables in 
order to keep the cost to a minimum. Thus it is soon that for 
a sinallor depth of the structure when the nonlinearity is more, 
thvj var3.ation in the pretension forces bot-ween the linear and 
nonlinear analyses is also more. As the depth increases, the I 
effect of nonlinearity decreases, and so does the change in 
the pretension forces. How-.^ver, there is no specific trend 
in the results obtained. 

6.4 fURTHBR WORK: 

An effort at the extension of the current work should 
be predominaxitly more practically oriented. The formulation 
for suoli a case should have at least two loading conditions, I 
a. ti.e dead load and the pretension in the cables, and b. the ; 
above loads plus live load. Depending on the live load chosen, | 
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various live load positions may "be necessary. The live load 
positions to give the worst live load effects may he found 
using the influence line diagrams. It nia:y he worthwhile to 
specify the thickness of the vfeh of box beams and to consider 
only the flange thickness as the design variable. The topolo- : 
gical parameters, especially the height of the tower and the 
points of support of the girders by the cables could also he 
taken as the design variables. ^ 

These additions would certainly make the formulation 
more comx^lex and the optimization procedure vrould he quite 
time consuming. The author feels that the Behaviour Model 
Method j76| as discussed in Section 4.3 would he an efficient 
approximation procedure to try and append to the synthesis 
ap'proach. Since multiple loading conditions arc involved, a 
iiord.iiioar analysis for each loading condition would be too 
time consuming. Alternately,- a certain reduction in the 
allov;ablo stresses may he made to take care of the nonlinearity; 
It is, however, felt that since tho optimum design arrived at 
through this effort would only be a skeletal streicturo , conside; 
ing only the linear elastic behaviour would ho justified. i 

An interesting point to observe from such studies is, j 

v/iajtiicr an optimality criterion like the fully constrained 
docign emerges out from these extensions also. In such a case,; 
a lormulntiori combining the resulting optimality criterion and I 
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the optimization procedure v/ith a reduced numlDer of variables 
would be very efficient. Such a procedure can also incorpora. te, 
with little effort, the availability of discrete thicknesses 
of plates, and the discrete variation of thickness of the 
flange within the length of a member. 

A study of the effect of the topological configuration 
on the cost of the structure would be very meaningful. The 
cables can be arranged in fan shape, harp shape or any other 
arbitra.ry shape. Study of the relative economy of these, 
similar to the one made by leonhardt and Zellner | 20j , would 
be very instructive. It would also be interesting to find the 
effect of the number of cables used on the optimum cost of the 
structure. Intuitively, it is felt tlmt the increase in the 
number of cables would decrease the total cost of materials 
used. . 

The synthesis approach and the program developed can be 
generalised and extended to the optimum design of guyed 
towers. This would req.uire the grouping of tho members having 
the same design variables and modification of the program to 
handle such a grouping. The three dimensional stiffness 
analysis procedure would have to be incorporated in the programi. 
The effect of large deformation should bo included- in the non- 
linear analysis procedure for this class of structures. The 
nonlinoarities in the behaviour of gu^red towers," both due to 
large doformation and due to beam- column effect are reported- 
to bo significant. 
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APPENDIX A 


SECANT STIFFNESS MATRIX AND ITS DERIYATIl/ES 

A.1 The Stiffness Matrix S for a "bending element in the 
element coordinate system is given as: 
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SEICC^+Cg) EI(C^+C2) 
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where and C^ for various cases are as follow: 
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l!l(C^+Cp) 


II n 

1 ^2 
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El(C^+Gp.) 


II r 
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(A.I) 


a) For zero axial force or v/hen axial force is ignored which 
leads to linear analysis 




= 4 


= 2 


(A. 2) 
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"b) For compressive axial force, 


0(Sin 0 — 0 Cos 0) 

0 

2(l-Gos 0) - 0 Sin 0 

(A. 5) 

0(0 - Sin 0) 

ri - - - 

(A. 4) 

2(l-Cos 0) - 0 Sin 0 

in which 0^ = v/here Q is the axial 

force in compression 


c) For tensile axial force, 
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in which 



cosh^ - sinhf) 
2(1-00311*) 4-^sinh1^ 

^ - y) _ 

2( 1-cosh^.) sinh¥ 


(A.5) 


(A. 6) 


3 ^ 

El 


where Q is the axial force in tension. 


For the purpose of this worh only cases a and h are 


req.uired. fhe fixed end moment under unifom load for an 
elemont in axial compression is given by, 


fixed end moment 


- 12 4 u^ tan u 


(A. 7) 


where, q, is the u.d.l. per unit length, and u - , 

A. 2 Computation of derivatives of the stiffness matrix with 

respect to the design variablesi 

Three typical expressions are involved in the elements of 
the stiffness matrix. The derivatives of these expressions are 
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computed, and the derivative of the stiffness matrix computed 
therefrom. 

T 4- T7l EA 

Let , Its derivative with respect to a 

variable is, 
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'A'lOV/ f 

0(Sin 0 - 0 Cos 0) .r. 

C. = — = ^ (A. 13) 

^ 2 ( 1-Cos 0) - 0 Sin 0 ^ 


SC^ 

‘5^* = 

How, 


D1 (Sin 0-0 Cos 0 + 0^ Sin 0) - Hi (Sin 0-0 Cos 0) 
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Assuming that the axial force Q does not change 
due to a change in I, 


M _ / 

E • 2i3/2 


21 


(A. 16 ) 


"61 obtained by substituting tho g_uantities from 

ii.l6 and A. 14 in A.12. 

Similarly, let, 
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(A. 17) 
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( J. C + SI Ids M 1 -il. 
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- (A.19) 
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Di(20 - 0 Cos 0 - Sin 0) - 172 (Sin 0.- 0 Cos 0) 
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(A. 20) 


6E- 

7?—* is obtained by substituting the q.uantities from A, 20 and A. 16 
in A. 18, 

and are calculated knowing the type of cross section 

and tho design variable. 



